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DECLARATION OF VISHWANATH R. IYER Ph D 
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I, VISHWANATH R. iyer, Ph.D., declare and state as 

follows : 

1. I am an Assistant Professor in the Section of 
Molecular Genetics and Microbiology, Institute of Cellular and 
Molecular Biology, University of Texas at Austin, where my 
laboratory currently studies global transcriptional control in 
yeast, gene expression programs during human cell 
proliferation, and genome-wide transcription factor targets in 
yeast and human. Immediately prior to this position. I spent 
four years as a postdoctoral fellow in the laboratory of 
Patrick O. Brown at Stanford University studying the 
transcriptional programs of yeast and of human cells My 
curriculum vitas is attached hereto as Exhibit A. 

2. Beginning in Dr. Brown's laboratory, where I 
helped to develop the first whole genome arrays for yeast and 
early versions of highly representative cDNA arrays for human 
cells, and continuing to the present day, i have used 
microarray-based gene expression analysis as a principal 
approach in much of my research. 

3. Representative publications describing this 
work include: 



scale. ••■scien« 27i;l80-6i"if^7°.;?"'? '""^nifi . 

using DNA inicroaT-ra,,o « « -"' ^"M- t-argec effects 
(1998) ^^'^''•^y^-v- .^Nature.Med..; 

the response of '^1'!;%^^''!^"^''^°"^^ in 
. . Science^S^Jfi^f^rrSsfy ^^^^ ---.^ 

expressr:„ variation in gene 

Two. Of , the papers describe our noo ^i: • 

^ our use of microarray-based 

(I=eR.si et al.. dur.ng.the shift from- fertnentation to - ■ 
respiration, .nd in hu^an cells ,l,er et al . ; hui^ ■ 
lioroblasts exposed to semmi ■ n 

use . ' • reference describes our 

and iTr""" '"^"'"^ '^-^ ' validation 

and .dentxf.cation of secondary dru. effects .Marton et ai " 

- one describes our use of expression profiling as a 

ec. discriminate ,a„on. bu„.n cancer 

cells (Ross et al.). ^aiicer 

4. Whether used to elucidate basic physiological 
responses, to study primarv ^r,H ^ o-^ogicai 

• . primary and secondary drug effects, or to 

discriminate and classifx/ >,„r„= 

Classify human cancers, expression profiling 
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Attached hereto as Exhibit B. 
Attached hereto as Exhibit c. 
Attached hereto as Exhibit D. 
Attached hereto as Exhibit E. 



as we have) practiced jt relies- fo^ • 

a££2~5 Of «^iessio«. - " r",:°"'=°^«i-" Of 

'5- -For example.,-, we have'deraonstrai-^H >w- 
use tSe presWibe-or abseS:ce^ of a I^K ' """""^^r^^"" «e can 

. • --f ■ ^ ® 9.^ a characteristic drnrr 

-signature- pattern of altered a.no ' 

cells to explore the .ech ion in drug-treated 

^.•^.or.e the mechanism .of drug .action anri ^« 
secondary, .effects that can sicm^V " ' i^'^ identify 

: ' ^ . . yr '. signal- poterit'iallV''a»=i««.-^- 

side effects. As an«r>,« ' ' ' - aeleterious drug 

AS another exainple, we have demonstrator! 
gene expression patterns t^an . ■ rated that 

. *'°^'-erns can be used to cls'sc-i f,. -\,- ■ 

cell .lines,, while it -i « > ' ■ - ' ^ "^^.^^^-."^^^^ tumor 

,Moxo.ica: funcj: oVL e:::::: •^-"'•^ - — - 

underlying these results tL "^-=har.i.sps, ^. 

retire .nowled.e or ^ 1 ■ '""""'^ - 

proteins ^olo.g.=,i function of, the. encoded 
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r: °' Pa^erns us.« in such 



comparisons is determined ^ the ^7Z:^::: 7T 
resolution of the pattern t*. 

resolution pettem^ ere^ene ali:" ^^'"^ 

■=b:^risons than loWer resoTu 

, . i^esolution patterns. with such v,,- k 

resolutions comes a pot-^^^ ^. higher 
stati.^• 1 .^°rrespondingly higher degree of 

statistical confidence for distinguishing different oatt 
as well as identifying similar ones. ^ 

7 . Each gene included as a probe on a mi n 
provides a signal that is specific to the 

at least to a first a.. • '° ^°^"^te transcript. 

rirst approximation * Esrh r,=, 

bach new gene-specific 

In a more nuanced view ir -i^ 

signal the presence of a variet" of sjl i ^riari-"'"'? ' ^"^^ ^<=> 

spiice variants of a single gene. 
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pbe,,,,dded to,a^ ^crcarray thus increases ' the' ni:;^".! g«,es 
detectable by the devicp •!t,^^«= • . aenes 

y ne device, increasing the resolving power of 

the device. : As I 'note.>abdve hioh*.!- ™ i ^. 

' '^^^her .resalution- patterns are • 
ge„eralix .„ore. useful ^^or^risons^^ acver resolution ' ' 
patterns., . Accordingly, each, new jene probe,..dded to a 
»icroarr.ay,:increase. .th, .usefulness of the device in .gene ' " ' 

expression P-rc«Ung .analyses,, .,,^is, proposition Is so wen- '' 

established as to be virtuallv . • ' ■ ■ . 

,.,oe yirt^ually . an jixiom i„^the,,art, and has 
been as.long as I, , have, .b,en forking in the field; -and' " ' 

oertainl^. since th, ti.e ,I en^ar.ed on the production 'o'f^ whole ^ 
gence arrays in early^l„,. singly put/ arrays witi, few^ - 
sene-specific probes are inferior to arrays with .ore gene- ' 
specific probes. , f 



e. fo-^ exa^le.. cur ability to subdivide , cancers 
into discriminable classes by expression nr„f ■ ' ^^"""^ 
by the resolution of the oat^ P-°*^ling .s limited 

ion Of the patterns produced, with more genes 
^ntributing to the expression patternsrwe can pot^Ju^- 
draw finer distinctions among the patterns, thus subdividing 
otherw.se indistinguishable cancers into a -greater number of 

likelihood that the cancer.! r-i ^ 

'=^^ssif led together will respond 
similarly to theraneutio -^espona 

y nerapeutic intervention, permitting better 
individualization of theranv =r,/i 

outcomes. "^"^ ' ^^"^^ ^"^^-"t 

9- If a gene does not change expression in an 
experiment, or if a oen^ •{ e 

a gene is not expressed and produces no 

(.Continued) 

without discriminating among them ar„q 

of a variety of allelic vaSa^tfif a sfniif ''"''^ '° "'^^^ '^^^ presence 
discriminating among them. single gene, again without 
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d.f^^^.^^?^^er:i9.e:^, is : not .to say "that " T"' 

r^er, of , conditions, have been' samnl.:,". ' 
... identify exDr*««<«« ^v. ^ ^ ?^?.A^^.r?an9>Aed to,:,,- 



P<f -^^^.7^ ®ivP^®ssion, changes . in fan*- . ' ' 

that a gene is no/' ^ ^ -^^^^^^^??'^./...an..experi^^ showing 

V, include — and sq a • -i »4i«uxa. 

->.... .... . ..,,- an,a fs. a. biologist. ., I would want mv 

i„_i„j^ . ■ '7''*f*^i'Pyiniicroarray to 

each newly ident If,- «^ °y co 

" ■ • ■ ■■ ■■ - • • ^-.^•^°?"5,i.f ied .gene as. a. prpbe. .• 



/ Ik • .t'. 



10. I declare further- Vv^^i:'' It i -^"-^^ ^ ^> 
^^rein- of cWn ^^l^l,^^ 
-de Oh ihfonnation and belief are b!l ! ^^atexnents 
further that these stat.. '° 
.that .Willful al e^ -de ,ith.,,the,3^o,l,dge. 

.punishable by; fine or 1 ^ '''^ '^ ^ 

. we or^ imprisonment, or both, under 

Section 1001 of Titl*» in 

"i. iicie 18 of .the United States oo,i« 
.jeopardize the validitv nf » states.^Cpde .and. may 

thi. a \ . ^"^^'^ °^ ?ny. patent application, in which 
this declaration is fii^^ ^ . -^*.:Wnica 

... n filed, or .any.patent that, is sues,,, thereon. 
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Vishwanath R. Iyer 

Assistant Professor 

Section of Molecular Genetics and Microbiology 
Institute of Cellular and Molecular Biology 

MBB3.212A, University ofTexas at Austin 

Austin, TX 78712-0159 

Phone: 512-232-7833 

Fax: 512-232-3432 

Email: vishy(S)mail.utexas.edu 

Education/Training 

Bombay University, Mumbai, India B Sc rioR-r-j rKor«;c^ c t>- i. • 

M. S. university of Baroda, Baroda. India u L^^i^^sS^J^'^''"""'^ 

Harvard University, Cambridge MA Ph.D. 1996) G^SS ^ 

Stanford University. Stanford CA Post-dortorii (^^^foooX Genomics 

Research Experience 

9/00-5/03 Assistant professor. Section of Molecular Genetics and 

Microbiology, University of Texas. Austin TX 

■ Global transcriptional control in yeast 

Sno^?'^?? programs during human cell proliferation 

- mZ™Hv ^^^"^^"Pt^°/ targets in yeast and human 
Loliaborative microarray facility 

5/96-8/00 Pojt-doctoral fellow Stanford University, Stanford CA 

(Advisor: Dr. Patrick 0. Brown) ^taniora LA 

■ Yeast whde-genome ORF and intergenic microarrays 
Human cDNA microarrays for expression profiling 

9/89-4/96 °3^--<|^-,3,^^,HaivardU„iversi.y. Cambridge MA 

■ Yeast transcriptional regulation 

Honours and Awards 

Government of India Biotechnology Fellowship (iQBy-iqRo^ 
University Grants Commission Junior Research Fellowship (1989) 
Stanford University/NHGRI Genome Training Grant (1996) 

Invited Conference talks (selected) 

Invited Lecturer, NEC-Princeton Lectures in Biophysics 
Pnnceton. NJ (June 1998) 

^911)'""°° 

'"^^^m-^^^^r" """'^""^ 



'°t?b&d (r^^^^^^^^^^ Conference 

Invited Speak»i^;^,,ucaora.Genor&(^iZ^I;^^'^ f . , 

Invited Speaker "BioArrays Europe:20Q3" o " ^003 

Cambridge, UK (Sep/Oct 2003) ' 




DepartmentaJ Seminars 

New York Uniyereity^choof of Medicine, ^ 
November 20 2002 

' ^^^r^"^"^ "^'"^ Seminar Series, 

UCW School of Medidne, I)epartment of Humaii Geietici " ^ 
June 2 2003 

National Human Genome Research Institute 

June 12 2003 " - i^v-^ 

Sangerlnstituteof the Wellcome Trust, Hinxton UK 
Sep 2003 ' 

Other Professional Activities 
Revi^ewer for Genome 

'"~'2?03 """^^ "^^^"^ using DNAMicroarrays" 

Member, NIDDK Special Emphasis Review Panel ZDKi (2001-2002) 

Publications 

^ -^'^'S' Poly(dA:dT), a ubiquitous promoter element that 
stimulates transcription via its intrinsic DNA stnicture, £mIo J.T4?257?.2S79. 

S t ?osl*7or«^ "'=3 TR and TC 



4.DeKslJ,L,:l3j6tii&Biw^P^'.^ . 



RoWts 



E. 0.& FriendS;,H. (1098) Dra£ taiS™M,h^ y-f'^"'?^ ^ 



6. Lutfiyya L L.,JvgrV. DeRisU „DeVit M J Brovm P n ii;'t«i;'^ 'l:>r"V 
Saccharomyces cerevisiae. Genetics 156:1377-1391 ^ regulate in 



8. l3jer\LR», Eisen M. B., Ross D. T SchulprP Ayr««^«T t , ^ ^ 
Staudt L. M., Hudson Jr. J BoLw M 1 ^li^? n" ct^ r 

Brown P. O. (1999) The iraAsSS^^ ^^'^^^ ^" ^^^stein D. & 

fibroblasts to seSaS283f^^^^^^ ^ 

9. De^i L & (.999) Genomics and array technology. Curr. Opin. Onool. 

10. Ross D. T., Scherf U., Eisen M B P*>rnn p m cr> n t> , 

Nature Genetics 24: 227-235 ^ expression patterns in human cancer cell lines. 

11. Sudarsanam P., lisiiV^ Brown p. O. & Winston F r2oon^^ 
expression analysis of sn/7su;i mutants nfc^f - ^ n ^-genome 

97: 3364-3369 Proc. Natl Acad. Sci .(USA) 



P-t«"'^^^^^^^ domain 
EMBOJ 19; 2323-2331 ATP-dependent chromatm-modifying factor 

13. Gross C, Kelleher M.. JyerV. R., Brown P. 0., & Winge D R f2ooo^ TH.ntifi ^ 
of the copper regulon in Saccharomyces cerevisiar^ DNAu^^r^^^^^^ 

Chem. 275: 32310-32316 inicroarrays. j. Biol. 

14. Reid J. L.. Tver V. R.^ Rrnwr. p q &ctnihlK^ roo^rv-* in j- ^ 



Kaiv^ J..:KoiiBo C WataiS^q Itph M.,-Shibata K. 

M., PeRisi J.L., Bken M^B Ue^^^^ Yc^hild A., Kusakabe . 



Immunoprecipitation on I^cSvs t^nM^^^"^^^ . 

Press, 2003). _ ' '^'^ 453-463 (Cold Spnng Harbor Laboratory 

*(not peer reviewed), \ , , ^ - ^ : . 

. : . . ■•■ . '.. : ' ■ ^ : ■ , \ i'. , ■ • 

19. Killion, p., Sherlock G. and hwV^ f2ooQ^T>,oT«^^ a 
open-source implementatio n of the W nrH ^iray Database, an 
Bioinformatics 4: 32 ^ ^ Microarray Database 

20. Hahn J. S., Hu Z., Thiele D. J. & iver V R r«««r« -iat- j . , 



Current/Pending Research Support 

NIH/NIAAA 

"INIA: Microarray Core" 

to define changes in gene exSon S ei heS^^^^^^ °' 

consumption. ""^^ P^^^'^^ or accompany excessive alcohol 

Role: Co-investi'gator 
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oi/bi/02 •^08/31/04 



J«?s Hi|hw Education Coordinating Board (ARPI 
._M.^y,teed global n«p^^^^^^^ ^ 

- -It 



This is apildt projcSrt to map the'ih vi^^ i -r ' 

Rcder-PI^^'^- • - 1.-, .,.„,.^^^..; ,. - ^. . ..n :v .; - 

NSF . - - . . ..„,>•• .■: ■ ■• •• I ■•■ • ■■ ■ 



"Feedback from Multi-Source Data Minincr r- • ^ . , 

Discoveiy" ^ - ^ . ^^.^ *° Expenmentation for Gene Network 

Role: Co-investigator ; . ' ' , , ; ^ ^ ^ ^ - ' 



iWo^- 'V30/T ^^^""^ "^'^ "'^"'^ 

NIH -■. 



"Analysis of genom^widetianscrip^^^ ' 

&DK.r^r ^ '"'""^'^^ ti^rough 
Role: PI 



1 ' 



us Amy Medical Research and Materiel Command 

hum» ceU. through the use o&X»Srrar&^ 

Texas Higher Education Coordinating Board f ATPl 

gSSSe"""^ hish-ft^ughput platform for measuring gene fi.nction on a 

h^^fSpufstJ^'c^as^ro^^lifxr ^ '^^-^ "'f »™ 

of gene funSon. "'^ ^"""""S "P'd and systematic evaluation 
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RKher-Vte. Sara 270. iB28(iggS). 

as. t; a Jam »ri a c. Bgia am cw a^. 3862 

•AO, CiB8Q:aPvovidO.S.Hogne&i.AocMstL>4eadi 
Stt USA aa. 263 (1991): .B. Ticfiench 

: CWBOJL 13. 3822 (1994); M.T.MB<fcK««,Crf 
' 87. 75 n 99ek D. a StohM. K. O. T«M. R. p. 

• '^AkttAaet Sot OSLAW. 7137 (1896). 



. 4068 0992): V. Sttvma. k. Suvwna. R Men- 
'* '^^.^ ^ HWfipetfi. J. escrm^ 174. 5057 

^ ^eaflarpt 178. 3015(1996), , 
37. U Ho cf a£. Caff 77. 869 (1994). 



UJetaI..09A)02UiK, filed Oe^^ 1997 

(PF-0459).;- . . =w# 



-Bchibit -B* attached to Dedaration of Vtshwanath 

R. Iyer. Ph.D. ' - ' * - ' 



Explwlnig thq Metabolic Snd Genetic^^^^ of 
Gene Expression on a Genomic Scale „ 

Joseph L DeRisi.^ 



^^^^^^f^''^^ every gene of Saccfta/omyces ce«Ws/ae were used 

to carry out a comprehensive investigation of the temooral nri,»L™!,7 • 
accompanying the metabofic shift from fSneStiK reSnn^Tl""''^'*'*'*'" 
p™fi.^ ot>served for genes wrth known me1Bctn"cJ°or;%SS"i^^^ 
metabolic reprogrammlng that occur durino the diauxic Shift i»nH fhlr ^^^^^^ °^ 
o^manypjeviouslyunchJacterizedgenesS^^^^ 

same DMA microarrays were also used to ^de^ti^ oenes who-sp «iX«ci!„ . « ".I ! 
by deietiori of the transcriptional co-repressor^S^^XSn of^^^^ 

nSfS"?" "'^■l; ^iemonstrate t^S Sl^rTnd X o^thT:? 

preach to genomewide exploration of gene .expression patterns. 



The complete sequences of nearly a doicn 
microbial genomes are known, arid in the 
next several yean we expect to loiow the 
complete genome sequences of several 
metaioans. including the . human genome. 
Defining the role of each gene in these 
genomes will be a formidable task, and un- 
derstanding how the genome fimctions as a 
whole in the complex natural history of a 
living organism-presents an even greater 
challenge.' 

Knowing when and where a gene is 
expressed often provides a strong due as to 
its biological role. Conversely, the pattern 
of genes expressed in a cell can provide 
detailed information about its state. Al- 
though regulation of protein abundance in 
a cell is by no means accomplished solely 
by regulation of mRNA, virtually all dif- 
ferences in cell type or state arc correlated 
with changes in the mRNA levels of many 
genes. This is fortuitous because the only 
sptciiic reagent required to measure the 
abundance of the mRNA for a specific 
gene is a cDNA sequence. DNA microar- 
rays, consisting of thousands of individual 
gene sequences printed in a high-density 
array on a glass microscope slide (/, 2), 
provide a practical and economical tool 
for studying gene expression on a very 
large scale (3-6). 

Sacchmomyces cerevisiae is an especially 
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favorable organism in which to conduct a 
^tcmatic investigation of gene expression. 
The genes are easy to recognize in the ge- 
nome sequence, cis regulatory elements are 
generally- compact and close to the tran- 
««iption units, much is already known 
about its genetic regulatory mechanisms, 
and a powerftil set of tools is available for its 
analysis. . . 

A recurring cycle in the natural history 
of yeast involves a shift from anaerobic 
^'f["C"^tion) to aerobic (respiration) mc- 
ubolism. Inoculation of yeast into a medi- 
um nch in sugar is followed by rapid growth 
tueled by fermentation, with the production 
of ethanol. When the fermentable sugar is 
exhausted, the yeast cells turn to ethanol as 
a carbon source for aerobic growth. This 
switch from anaerobic growth to aerobic 
respiration upon depletion of glucose, re- 
ferred to as the diauxic shift, is correlated 
with widespread changes in the expression 
of genes involved in fundamental cellular 
processes such as carbon metabolism, pro- 
.i^"^"'** and carbohydrate storage 
(/;. We used DNA microarrays to charac 
tcme the changes in gene expression that 
take place during this process for nearly the 
entire genome, and to investigate the gc- 
iietic circuitry that regulates and executes 
this program. 

Yeast open reading frames (ORFs) were 
/Snox polymerase chain reaction 

(PGR), with a commercially available set of 
primer pairs (8). DNA microarrays. con- 
taming approximately 6400 distinct DNA 
sequences, were printed onto glass slides by 
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-iwing a simple robotic printing device (9). 
Cclb from an exponentially growing culture 
of yeast were inoculated into ftwh medium 
«nd;&rown at 30^C for Zl houti. .Aftcr;an 
initial 9 hours of growth, samples were har- 
vcsted at^sevcn successive 2-hour intervab, 
and mRNA was isolated ( JO). Fluorcscently 
labeled cDNA was prepared by.revcise tran- 
scnption in the presence of ey3(grecn)- 
or 'Cy5(red)-Ubclcd dcoxyuridine triphos- 
phate (dUTP) (J /) and then hybridized to 
the microarrays (12). To maximize Ae re- 
'liability with which changes in expression 
levels could be discerned, we labeled cDNA 
prepared- from celb at each successive time 
point with Cy5. then mixed it with a Cy3- 
labeled "reference" cDNA sample prepared 
from cells han^ested at the fim interval 
after inoculation: In this experimental de- 
sign, the relative fluorescence intensity 
measured for the Cy3 and Cy5 floors at 
each , array clement provides a reliable mea- 
sure of the relative abundance of the corre- 
sponding mRNA in the two cell popula- 
tions (Fig. I). Data from, the series of seven 
f consisting of more than 
43,000 expression-ratio measurements, 
were organized into a database to facilitate 
efficient exploration and analysis of the 
results. This darabase is publicly available 
on the Internet (13). 

During exponential growth in glucose- 
rich medium, die global pattern of gene 
expression was remarkably stable. Indeed, 
when gene expression patterns between the 
first two cell samples (harvested at a 2-hour 
interval) were compared. mRNA levels dif- 
fered by a factor of 2 or more for only 19 
genes (0.3%). and the largest of these dif- 
ferences was only 2.7-fold ( M). However, as 
glucose was progressively depleted from the 
growth media during the course of the ex- 
pcriment, a marked change was seen in the 
global pattern of gene expression. mRNA 
levels for approximately 710 genes were 
induced by a factor of at least 2, and the 
mRNA levels for approximately 1030 genes 
declined by a factor of at least 2. Messenger 
RNA levels for 183 genes increased by a 
factor of at least 4, and mRNA levels for 
203 genes diminished by a factor of at least 
4. About half of these differentially ex. 
pressed genes have no currently recognized 
function and are not yet named. Indeed 
more than 400 of the differentially ex- 
pressed genes have no apparent homobgy 



to any gm whose function is known (/5). 
The responses of these previously unchar- 
acterittd geries to the d'tauxic shift therefore 
provides the first small cliie to their possible 

The global view of changes! in expres- 
sion of genes with knowiv functions pro- 
vides a Vivid picture of the way in which 
the cell adapts to a changing . environ- 
ment. Figure 3«hows a ponion of the yeast 
metabolic pathways involved in carbon 
and enerjgy metabolism. Mappirig the 
changes we observed in the mRNAs en- 
coding each enzyme onto this framework 
allowed us to infer the redirection in the 
flow of metabolites through this system. 
Wc observed large iriductions of the genes 
coding for the eniymes aldehyde dehydro- 
genase iALD2). and .'acetyl-coemymc 
A(CoA) synthase (ACS J), which func- 
tion together to convert the products of 
alcohol dehydrogenase into aceryl-CoA, 
which in turn is used to ifiiel the tricarbox- 
ylic acid (TCA) cycle and the glyoxylate 
cycle. The concomitant shutdown of tran- 
scription of the genes encoding pynjvatc 
decarboxylase and induction of pyruvate 
carboxylase rechannels pyruvate away 
from acetaldchydc, and instead to oxalac- 
etatc, where it can serve to supply the 
TCA cycle and gluconeogcnesis. Induc- 
tion of the pivotal genes PCJCi. encoding 
phosphocnolpyruvate carboxykinase, and 
FBPl, encoding fructose 1,6-biphos- 
phatase. switches the directions of two key 
irreversible steps in glycolysis, reversing 
the flow of metabolites along the revers- 
ible steps of the glycolytic pathway toward 
the essential biosynthetic precursor, glu- 
cosc-6-phosphatc. Induction of the genes 
coding for the trehalose synthase and gly- 
cogen synthase complexes promotes chan- 
neling of giucose-6-phosphate into these 
carbohydrate storage pathways. 

Just as the changes in expression of 
genes encoding pivotal cmymes can pro- 
vide insight into metabolic reprogram- 
ming, the behavior of large groups of func- 
tionally related genes can provide a broad 
view of the systematic way in which the 
yeast cell adapts to a changing environ- 
ment (Fig. 4). Several classes of genes, 
such as cytochrome c-related genes and 
those involved in the TCA/glyoxylate cy- 
cle and carbohydrate storage, were coord i- 
nately induced by glucose exhaustion. In 
contrast, genes devoted to protein synthe- 
sis, including ribosomal proteins, tRNA 
synthetases, and translation, elongation, 
and initiation factors, exhibited a coordi- 
nated decrease in expression. More than 
95% of ribosomal genes showed at least 
twofold decreases in expression during the 
diauxic shift (Fig. 4) (13). A noteworthy 
and illuminating exception was that the 



genes encoding mitochondrial ribosomal 
genes were generally induced rather than 
r^ressed. after glucose limiution. high- 
halting the requirement for mitchondrial 
- biogenesis ( 13). As more, is learned about 

every* gcrie in the yeast 
genome, the.ability to gain insight into a 
cell s response to a changing environment 
through. its global gene expression patterns 
will become tncreasirigly powerful.- ; 

Several distina temporal patterns of ex- 
pression could be recognized, and sets of 
genes, could be .grouped on the basis of the 
similantics in their expression panems. The 
characterhed; membcris of each of these 
groups also shared, important similarities in 
their hmaions. Moreover,, in most cases, 
common regiilatory mechanisms could be 
infened for sets of genes with similar expres- 
sion profiles. For -example, seven genes 
showed a late induction profile, with mRNA 
levels increasing by more than ninefold at 



Reports 



dw last^imqxjint but less than duc^U at 
thnepbiitt jl^g. -5B). AD of 
were :kribwn: to, be giucoK*rc- 
preoed, and five of the seven' were previously 
, '^^'ted to ^lare a comnxm t^stii^^ 
inig scquerice nJAS), the carbon source re- 
spmse elcmdit (CSREIX 16-20).. A search 
in the promoter regions of the remaining two 
genes, ACK/ and IDPl repealed that 
ACf?i, a ^e essential fior ACSl activity, 
also possessed a consexuus CSRE motif but 
interestingly. /DP2 did not. A search of Ae 
crinre yeast, gervome sequence for the con- 
scrisus.CSRE motif revealed caily four, addi- 
tional candidate genes, none of which 
showed a similar' induction; 

Examples ' from additional groups of 

expression profiles are 
illustrated in Fig. 5, C through F.. The 
sequences upstream of the named genes in 
Fig: '5C all contain stress response -ele- 
ments (STRE), and with thd exception 



fi 




Rq. 1. Yeast genome microarray. The actual size of the microarra!^n^!IIT^r^ ^ 
nwcroartay was printed as described roi tw. "iiwoarray is 18 mm by 18 mm. The 

scanning 'confocS mi^o^oS u^o ?ollI5 a UhfdaTJ"" ''^m^ "™ 
CONA probe was prepared tr^ m1^I°isSS U^J^L^e^^^Vn^'T*" ^ 

<0.2 g/lrte,) by reverse transcriptionln th™?^^^^^^^ ^"I" ^'"^^ ^ 

CyS-dUTP-labeted cDna nhai k mBwi li!; . Cy5-dUTP. in this image, hybridization ot the 
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of HSP42, have previously been shown to 
be controlled at least in part by these 
elements {21^24). Inspection of the se. 
quences upstream of HSP42 and the two 
uncharacterized genes shown in Fig. 5C, 
YKL026c, a hypothetical protein with 
similarity to glutathione peroxidase, and 
YGRCM3c, a putative transaldolase, re- 
vealed that each of these genes ako pos- 
sess repeated upstream copies of the stress- 
responsive COCXrr motif. Of the 13 ad- 
ditional genes in the yeast genome that 
shared this expression profile [including 
HSP30, ALD2. OM45, and 10 uncharac- 
teriied ORFs (25)]. nine contained one or 
more recognizable STRE sites in their up- 
stream regions. 

The hecerotrimeric transcriptional acti- 
vator complex HAP2,3,4 has been shown 
to be responsible for induction of several 
genes important for respiration (26^28). 
This complex binds a degenerate consensus 
sequence known as the CCAAT box (26). 
Computer analysis, using the consensus se- 
quence TNRrrGGB (29), has suggested 
that a large number of genes involved in 
respiration may be specific targets of 
HAP2,3,4 (30). Indeed, a putative 
HAP2,3,4 binding site could be fourd in 
the sequences upstream of each of the seven 
cytochrome c-related genes that showed 
the greatest magnitude of induction (Fig. 
5D). Of 12 additional cytochrome c-relat^ 
genes that were induced, HAP2,3,4 binding 
sites were present in all but one. Signifi- 
cantly, we found that transcription of 
HAP4 itself was induced nearly ninefold 
concomitant with the diauxic shift. 

Control of ribosomal protein biogenesis 
is mainly exerted at the transcriptional 
level, through the presence of a common 
upstream-activating clement (UAS ) 
that is recognized by the Rapl DNA-bin3- 
ing protein (31, 32). The expression pro- 
files of seven ribosomal proteins are shown 
in Fig. 5F. A search of the sequences 
upstream of all seven genes revealed con- 
sensus Rapl -binding motifs (33). It has 
been suggested that declining Rapl levels 
in the cell during starvation may be re- 
sponsible for the decline in ribosomal pro- 
tein gene expression (34). Indeed, we ob- 
served that the abundance of RAP J 
mRNA diminished by 4.4-fold, at about 
the time of glucose exhaustion. 

Of the 149 genes that encode known or 
putative transcription facton. only two, 
HAP4 and S1P4 , were induced by a factor of 
more than threefold at the diauxic shift. 
SIP4 encodes a DNA-binding transcrip- 
tional activator that has been shown to 
interact with Snfl. the "master regulator" of 
glucose repression (35). The eightfold in- 
duction of S/P4 upon depletion of glucose 
strongly suggests a role in the induction of 



downstream genes at the diauxic shift. 

Although most of the transcriptional 
responses that we observed were not pre- 
viously known, the responses of many 
genes during the diauxic shift have been 
described. Comparison of the results we 
obtained by DNA microarray hybridiza- 
tion with previously reponed results there- 
fore provided a strong test of the sensitiv- 
ity and accuracy of this approach. The 
expression patterns we observed for previ- 
ously characterized genes showed almost 
perfect concordance with previously pub- 
lished results (36). Moreover, the differ- 
ential expression measurements obtained 
by DNA microarray hybridiiation were re- 
producible in duplicate experiments. For 
example, the remarkable changes in gene 
expression between cells harvested imme- 
diately after inoculation and immediately 
after the diauxic shift (the first and sixth 
intervals in this time series) were mea- 
sured in duplicate, independent DNA mi- 
croarray hybridizations. The correlation 
coefficient for two complete sets of expres- 
sion ratio measurements was 0.87. and for 
more than 95% of the genes, the expres- 



sion ranos measured in diese duplicate 
cxpenments differed by less than a factor 
of 2. However, in a few cases, there were 
discrepancies between our results and pre- 
vious results, pointing to technical limita- 
tions that will need to be addressed as 
/ir ^^^o'ogy advances 

u k L ^^"^ exceptions, 
the high concordance between the results 
wc obtained in these experiments and 
those of previous studies provides confi- 
dence in the reliability and thoroughness 
of the survey. 

The changes in gene expression during 
this diauxic shift are complex and involve 

L®'^ °^ information 

about the nutritional and metabolic state 
of the cell. The large number of genes 
whose expression is altered and the diver- 
sity of temporal expression profiles ob- 
served m this experiment highlight the 
chal enge of understanding the underlying 
regulatory mechanisms. One approach to 
defining the contributions of individual 
regulatory genes to a complex program of 
this kind is to use DNA microarrays to 
Identify genes whose expression is affected 



Rg. 2. The section of the ar- 
ray indicated tyy the gray box 
in Fig. 1 is shown tor each of 
the experiments described 
here. Representative genes 
are labeied. In each of the ar- 
rays used to analyze gene 
expression during the diauxic 
shift, red spots represent 
9enes that were induced rel- 
ative to the initial timepoint. 
and green spots represent 
96nes thai were repressed 
relative to the initial timepoint. 
In the arrays used to anafyze 
the effects of the tuplA mu- 
tation and YAPi overexpres- 
sion, red spots represent 
genes whose expression was 
increased, and green spots 
repnesertt genes whose ex- 
pression was decreased by 
the genetic modification. Note 
that distinct sets of genes are 
induced and repressed in the 
different experiments. The 
cornptete images of each of 
these arrays can be vtewed on 
the Intemel (73). Cefl density 
as measi^ed t)y optical densi- 
ty (OD) at 600 nm was used to 
measure the growth of the 
culture. 
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that mult ftom deletion of die TLfPi gene nmaUmofZt^\^i^ '^^^ 
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glucose requires die DNA-binding repressor (40), 
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R9' 3. Metabolic reprDgrammrno inferred from alobai anaiueie ^« 

•netatxfc irterrediat^aBBeSlfiedX^^^ ^P^'""- Only key 

in this metabolic ciicult are identified tv ram^ttS^^T^S^^ "'^^ ««P 

and flIycogefHWjranching enzyrne SwSSi^E^^ 
VPR1 W show significant hom^ toJm«!m rtentdied. but the ORFs YGR244 and 

White letterine ide^Wy T^I^^^S^^^^^^^^^^^ ^ 
green lelterinflSy'genes whose^^L^^X 5f 

induction or repression is indicated ^^V^^^^^^J^^^'^ ^ "«9nitude of 
succinate dehydrogenase, the indicated ^rSS^ uS^""'**"'' 

genes listed in the box. Black and white boxes in^no ^^^^ni^^^^** °' ^" 

twofold). The direction of the arn>ws connectSSwT^S. h'*^' ^^'^ <^ ««" 
flowof metabolic ntermediates. inferred from tS^^^L^^J^t'^^r'^'^*^ *'''«^'p" 
rapresenthg steps catalyzed by Sen^^lS^^^^t^J^J^ ''l!^. 
The broad grey arrows represent major increa^TSlZ/L?^J2^f^ ^ h-ghlghted in red. 
inferred from the indicated changes h^ne^^^^ ^ "^^^ a««er the diauxic shift. 
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Wild-type yeast celb and cells bearing 
^tipn qf the^TUPJ gene (mplAwf 
grown in parallel cultures in rich medium 
contaming glucose as the carbon source. 
Messenger. RNA was isolated frohTeS 
nentialiy growing cells from the two poo- 
ulations and used to prepare cDNA la- 
beled with Cy3 (green) and Cy5 (red) 
respectively (i i). The labeled probes were 
mixed and simuluneously. hybridiied to 
the microarray. Red spots on the mictoar- 
lay therefore represented genes whose 
transcription was induced in .the.. tid>] A 

presumably repressed by 
i upl Mi ). A representative section of the 
microanay _(Fig.. 2, bottom middle panel) 
Illustrates that the genes whose expression 
was affected by theiupiA mutation, were, 
in general, distinct from those induced 
upon glucose exhaustion (complete images 
of all the ansys shown in Fig. 2 MK unA. 

34 f l^W k'"" Neverthelea. 
34 (10%) of the genes that were induced 
by a factor of at least 2 after the diauxic 
shift were similarly induced by deletion of 
TUFl, suggesting that these genes may.be 
subject » TUPJ-mediated repression by 
gUjcose. For example, SUC2. the gene ^ 
coding invertase. and all.five hexose trans- 
porter genes that were induced during the 
course of the diauxic shift were siinilarly 
induced, .n duplicate experimena. by the 
deletion of TUPi. , j 

The set of genes affected by Tupl! in this 
wcperiment also included a-glucosidases. 
the.mating-type^pecific genes MFAJ. and 

RNR2 ar^d RNR4, as well as genes involved 
in flocculation and many genes of unknown 
function. The hybridization signal cone- 
spending to expression of TUPl itself was 
also severely reduced because of th^ (in- 
complete) deletion of the transcription unit 
m the tup/A strain, providing a positive 
control in the experiment (42). 

T..I^'?M°^ transcriptional targets of 
Tupl fell imo sets of genes with related 
biochemical functions. For instance, al- 
though only about 3% of all yeast genes 
appeared to be TUPl -repressed by a factor 
of more than 2 in duplicate experimena 
under thew conditions. 6 of the 13 genes 
that have been implicated in flocculation 
[15) showed a reproducible increase in 
expression of at least twofold when TUPl 
was deleted. Another group of related 
genes that appeared to be subject to TUPJ 
repression encodes the serine-rich cell 
wal mannoproteins. such as Tipl and 
Iirl/Srpl which are Induced by cold 
shock and other stresses (43). and simiUr. 

the seripauperins 
m. Messenger RNA leveU for 23 of the 
^6 genes in this group were reproducibly 
elevated by at least 2.5-fold in the nipiA 



....... 

nrrain ..Jio-r.i ... >->iM. -..-^.ir...^.. 



"ram, and 18 of chew ^nes were induced 
by more than «eyenfold when TUPl was 
cqnaait.' none of 83 genei that 
opuM be cla»ified a« putativeregulatora of 
the cell divition cycle were induced more 
than , twofold by deletion of TUP; Thus 
despite the diversity of the. re^ktoty sys-' 
terns that employ .Tiipl , most of the geiies 
that It reguUti» under these, ^itions 
rail into a h'mited number of distinct func- 
tional classes: . ■, •, 

^Because the microarray allows , us to 
moniwr expression of nearly every gene in 
yeast, we can, in principle,, use this ap- 
proach K> identifyrall the. triucriptional 
targets of a regulaiory.protein like Tupl. It 
IS important to note, hpwever, that in any 
single experiment of this kind-we can only 
recogniie those urget genes that are nor- 
mally repressed- (or induced) under the" 
conditions of the experiment. For in 
stance, the experiment described here an- 

and MFA2, the genes encoding the a- 
taaor mating pherorobne precursor, are 
normally repressed. In the isogenic tubl^ 
stnim. these genes' were inappropriately 
expressed, reflecting the role that Tupl 
plays in their repression. Had we instead 
canted out ,a»ii experiment with a MATA 
strain (in which expression of MFAl and 
MFA2 IS not repressed), it would not have 
been possible to conclude anything re- 
siding the role of Tupl in the repression 
of these genes. Conversely, we cannot dis- 
tinguish indirect effects of the chronic 
ateence of Tupl in the mutant strain from 
effects directly attributable to its panici- 
pation in repressing the transcription of a 
gene. 

Anodier simple route to modulating the 
activity of a regulatory factor is to ovetex- 
Pt«s the gene that encodes it. YAPl en- 
codes a DNA-binding transcription factor 
belonging to die b-2ip class of DNA-bind- 
ing proteins. Overexpression of YAPJ in 
yeast confers increased resistance to hydro- 
gen peroxide, o-phenanthroline, heavy 
metals, and osmotic stress (45). We ana- 
lyied differential gene expression benveen a 
wild-type strain bearing a control plasmid 
and a strain. wid» a plasmid expressing YAPJ 
under the control of the strong CALLIO 
promoter, bodi grown in galactose (that is 
a condition that induces YAPJ oveiexpres-' 
sion). Complementary DNA from the con- 
trol and YAPi overexpressing strains, la- 
beled wiA Cy3 and Cy5. respectively, was 
prepared from mRNA isolated from the nvo 
strains and hybridised to the microarray. 
Thus, red spots on die anay represent genes 
iM V/Si " overexpress- 

Of the 1 7 genes whose mRNA levels 
increased by more than dweefold when 



2f .i was overexpressed in this wty, five 
bear homology to aryl-alcbhol oxidoredue^ 
tases (Fig. 2.and Table l):i An -addi^ 
four.ofd>e.ger|es-in this set also beloiab 
Ae .general dass of dehydrogenaseSti- 
doreductaws. Very little is fawTbout 

""doreductases in - 
5.,cerew««.,but these'eniymes have. been 

hev '^e*' '-hich. 

tionsi. oxiduing aromatic, and aliphatic 
unsaturated alcohols to aldehydes wid, the 
production of hydrogen peroxide (46. 47) 
TT>e,fac,_that a remarkable fraction trf^ 
"tgcts.identifiedKin.this experiment 
long to the same small, ftinaional group of: 
oxidoreductases suggests that die^ g^es 
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r w uiese genes encode surar bcr. 

We sarehed^ j[!P'-*''«J»K «tcs 
(TTACn^AA or TGA(n^AA) in the «^ 
quences^upmeam of the laiget gene, we 

induced by more tha^ 

i^^ ^T ^^Pi ovenaUsiohS 
««J* more.bi„ding rites widiin 600 b^' 

of the start codon (Table 1). 
i«mg that Aey are direcdy reguUtJ^ 
rapl. The absence of canonical YapUbind- 
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R9.'4:'CoorrfinatecJfB9- •* 
■ ulation of ;functionally re- 
■ lated genes. The curves 
represent the average irv 
duction a repression ra- 
tios for all the genes in 
each indicated group. 

total rttjmber of 
genes in each group was 
as follows: ribosoma) 

proteins. 112; translation 

etongatiofi and initiation - 
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reductase 

Transcnptionai activator involved in 
oxidative stress response 

Homology to aryl-aJcohoI 
dehydrogenases 

Putative glutathione transferase 
"Kpt^ ctehydrogenase 

Putative aryj-aicohol reductase 
Ammotriazote and 4-nitroquinoline 
resistance protein 

Homology to benomyl/methotrexate 

resistance protein 
Hypothetical protein 

NAPDH dehydrogenase (old yellow 
enzyme), isoform 3 

" v?2^'° 'VPOthetical proteins 

YCRl02candYNL134c 
Homrtoar to hypotheiical protein 

tMR251w 

NAD{P)H oxidoreductase (old yellow 
enzyme), isoform i ^ 

Similarity to A. thaliana zeta^rystaDin 

nomolog 
Malate dehydrogenase 
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Fold- 
icreas 

12.9 
10.4 
9.6 

9.0 

7.4 
7.0 

6.5 
6.5 

6.1 

6.1 
6.0 
5.8 

4.7 

4.5 

4.1 

3.7 

3.3 



ing,5itcs iqxmam of the othen may rtflea 
m ability of Yapl to, bind «itci that differ 
ftom the caiMmical biiding sites, perhaps in 
cooperation widi odier factors, or lew like- 
iFt n?ay rqjiescht aii indirect effect of Yapl 

one -or' more 

uitenTOdiary ^ctbn. Yapl sites were found 
otHy four times in the coireaponding region 
of an arbitrary tet of 30 genes that were not 
ditterentially r^lated.by Yapl.; , ' ; 

, Use of a D)>IA 'microarray to characteN 
ue the transcriptional consequences of 
muutions affecting the activity of reguia- 
tory molecules provides a.simple and'pow- 
eifiil approach to dissection and character- 
ization of regulatory pathwayis . and net- 



worb,.This saategy aUo.has an important 
practical application in drug- screening. 
5*utauons in specific genes encoding can- 
didate dnig urgets tan wrve as surrogates 
for the^ideal chemical inhibitor or modu- 
fawr pf.their activity. DNA.microanays 
can be .used to define the resulting sigha- 

ture pattern of aherations .in gene expres- 
«pn. and then subsequently used in an 
assay to screen for compounds -that repro- 
^""^ signature pattern. 

UNA microanays provide a simple and 
economical way to explore gene, expres- 
sion patterns on a genomic , scale. JThe 
hHrdles to exrending, this approach to any 
other.prganism are minor. The equipment 
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Rg. 5. Kstind temporal patterns Of ifxJuction or fBDrBS!i« twi« 

Poverties. (A) Temporal proNe of the cell denST?^J^^°"^*"** ^^'^ '^"^'^^ 
coocemratlof^ in the media. (B) Sevan genes ^uS^fJ^!^ OD at 600 nm and glucose 
the last tlmepcmt (20.5 hows). With the^SfeS^S^^r^li^ "^"^ °^ at 

were r» addrtional genes observed to match ttepmteffi^^ has a CSRE UAS. There 

by earty induction with a peak in mRNA levels at^hSf^T!!^"' of a dass of genes mart<ed 
repeats in their upstream promoter regions. (0) Cviodv^TV? S ® STRE motif 

reductasegenes. Marked byanlndu«S?SnSeSS^ c 
a consensus binding motif lor the HAP2.3 rnrn/lTcnmfi^^ 

expression prolile. (E) S4MI. GPPl. and s««^SL^^' shared a similar 

diauxic shift, and comhue to be repressedurcn^ V..^^ ^pressed before the 

genes comprise a large dass of genes that ar^m.^^ T^^T' Protein 

profited here contains one or mcie^i -bK^^.^^ f °' ^<=^ °' Senes 

tional regulator ol most rixtsomal proteins^ * "^'^ °' promoter. RAPi is a iraracrip- 



9 



www.sciencemag.org • SCIENCE • VOL. 278 . 2^ OCTOBER 1997 



•wjuired for Abricating and usiiis ;mA 
tn'croarray, (9) consim: of comJoJ^ 

.mphcty. It •as feasible for . small 

. 6«» genes to aboui 4 mohthi and 
once the amplified gene sequences were in' 
^ only 2 days, were required to print a 
set of no ratCToarrays of 6400 elerhents 
eacft.^ rtobe preparation, hybridiiitiori- 
and. fluoresant imaging are also simple 
procedures. Even conceptually simple oi- 
periments. as we described here, can yield 
^t amounts of information. The valdc of 
the information f^om each experiment of 
this kind will progressively increase as 
more « learned about the functions of 
each gene, and as additional experiments 
define the global changes in gen^expres- 
sion rn divene other naniral processes and 
genetic penurbations. Perhaps the greatest^ 

m/rK^^ " ?'*«'°P ■ «ff«cient. 
methods for organizing, distributing, iriter- 

pretmg. and extracting insights from the 
te'r °' 

. ' ' f ■ 

REFERENCES AND NOTES 

3. 0. Lashkart.ftoe. Atett Aeaa Sd U&Jk. kii,^ 

4. J. Defiisi «ten« <^ lX45V%S)'^ 
(X^**** « -t. /V«un, ao*^^ ,675 

f ■ li" * "■- Science 274, 610 (1996). 

Spnr« Hvtor Laboratory Press. (ST SoS» 
bor.r«'.l992),p.l93. «>pnnB Hw- 

ft^>«to each knom, or predicua proWn codta 
S^l^a wim the proioool swpS^ 

sw^MC as a tampiaie Each pw>ua\S 
«»^by agarose get el8ctrophor«rSo 

oeOTM cofteet It the lane contained a^rS.5S 

-n^ovwaOsuccesswelSB^^ 
pass ampHcaten of 61 1 6 ORFs was -94 JSi^ 

Glass sMes (Gold Se* were oeaneo for 2 toui to . 
soMion Of 2 N NaOH and roSeSSd^iS^^ 

tor 1 hoi,, and men dhad lor 5 n*, UW^^ 
^ouuTi o««a Orw samples liwn lOO-ol PC»»I^^ 

etnanol purifcatiSn^rBe!^^ 
"~«*terpa«. T)* resulting precipitates «^ 
««Pended m 3x standard sttine titrate <ssei «3I 

^^t^^^ *° » batd^oi?^ 

*°|«"«a^Oedu/pbrown. After prM. 
no^mooaroys were rehydraied for 30 s ina 
tMWd Chants and ttH« snap^ tor 2 ion 1^,2 
Wte noO-Q. The DNA was ^M^wt 

»«sW|adtoihesurtaceby8ubjeelt,8^S1i^^ 

polyH-lysine surface was Uoeked bv a i&^ 
";««Mtion in a sowion of 70 ntf^ suedr* ar<>^ 

n^2-pytroWinone (AkMch) 35 rrtol 1 M 
bone aod <pH 8.0. CJirectly alter the btocWng rie- 



e 



685 



• lion, thebound ONA wm dsrwtmdtyy 

. qjbatkm n (SsOed vvater St ^ T>« sides wve 
, vwn transfonsd hto s bsth of 1 00% etrmi St fcm 
■^tampentue. rinnd. snd then spui fri • cMctf 

uOerMiQi Ste wm ttofod in s dCMd bot« 

• nxvn tarnperatm irtl used. 

la .YPO medkim (8 itsrs). In s. 104ler te niieii u nk»i 
vessel, inbcUsted vvfth 2 ml of'a'i^ over- 

nigm cutu0 ct yeast strain DBV7286 (MATa. ua3, 
GAL2}. The fermentor was rnaintahed at 30*C wmi 
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from the fermentation vessel by peristaitic purr^. 
spun down at room temperature, and then flash 

trozen with liquid nltrogea Frozen cells were stored 
at-80X. - 

rated duhng reverse transcription of Viis iig of 
po^^denytated |poly(A)*I RNA. primed by a 0106) 
«*9omer. This mixture was heated to 70*C tor 10 " 
min. and ttwi transferred to ice. A premixed solu- 
tiorv, consisting of 200 U Superscript II (Gibco), 
^er, deoxyribonucleoside triphosphates, and flu- 
orescent nucleotides, was added to the RNA. Nu- 

deotides were used at these final concentrations' 
500 for dATP, dCTP, and dGTP and 200 
lor dTTP. Cy3-dirP and CySnSL/TP were used at 
a final concentration of 100 »tM. The reaction was 
then incubated at 42*C for 2 hours. Unincorporat- 
ed fjuorescem nucteotides were removed by ftrst 
dilutirtg the reaction mbrture with of 470 »d of 10- 
mM tris-HO (pH 8.0)/i mM EDTA and then subse- 
quently concentrating the mix to *5 jiJ, using Cer>- 
tricon-30 mi croconceniiators (Amiccn). 

12. Purified, labeled cONA was resuspended in n,Jcxf 
3.5X SSC containing 10 tig polyfdA) and 0.3 »j of 
10% SOS. Before hybridization, the solution was 
boied for 2 min arxJ then allowed to coot to room 
temperature. The solution v^ applied to the mi- 
croarray under a cover slip, and the slice was 
placed in a custom hybridization chamber wttch 
was subsequently incubated for -8 to 12 hours in 

a water bath at 62*C. Before scanning, sfides were 
v«shedin2xSSC.0.2%SDSfor5rrtn.anothBn ' 
0.05X SSC for 1 min. Sfides were dried before 
scanning by oentrifugation at 500 rpm n a Beck- 
man CS-6R centrifuge. 

13. The comptete data set is avatebte on the Internet at 
cmgmstantord.edu/lpbrowrVexplore/index.html 

14. For 95% rtal the genes analyzed, the mRNA levels 

Pleasured in cells hanested at the first m second 
interval after irK)cUation diftered by a factor of less 
than 1 .5. The conel aii on coefficient for the compff- 

ison between mRNA levels measured lor each gene 
in these two dmererrt mRNA samples was 0.98. 

VSffwn Ajpficate mRNA preperations from the same 
ctfsampte were compared in the same way, the 
correlation coefficient between the expressnn levels 
measured lor the two samples by comparative hy- 
txidization was 0.99. 

15. The numbers and identltjes of knom and putative 

genes, and their homologies to other geries. were 

gathered from the totowing public databases: Sac* 

charori^cBS Genome Database tgenome-www. 

stantord-ecM. Veast Protein Database (quest7; 

proteome^corn). and Munich Information Centre for 

Protein Sequences (speedyjnipsi)iochem/rv)g.de/ 

mips^east/index.htrnb}. 
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H^JS^^ P«*e. was defin«j.as haveig.pn 

^^25? ^ ^ ^^>^ « ^fi^ howa «id 
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2a-& LForstxjgandU Guarente.-Gar»es iW. 3,ii€6 
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29. Sngle^er abbreviations tor the amno add resi- 

^ ^* ^' ^* Asp; E cm: F. 

^'^ ^ ^y^' U Leu: M. Met N. 
. ?S^iL^V°',2!f '^'^-^SenT,Th^ 

B-C. G. or T; W-G. A. T. or C; R^ or G;.8nd YrC or 
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"'?°* t^?^!^ ^ Kalogerxipai^, Conw Appl. 
..^ Bosa 1Z 363 (1996). 

?i* d' ^""^ <06 (1994) 

^ iwd 4. 64 (1988). 
' ' ^^oegeneiate consensus sequence WCYRNNC- 
MW «^ used to search tor potential RAPi .btx*w 

^ as aioweu. 

Cee,ScL 15.3187(1995). 

, ;36.. For example, we observed targe inductions of the 
genescodjjg fof PCK7. FBPi (z. et a/.. Mo/. 

^ ^^^ W^tJ'^ "^^^ Styaxytate 
cy&e gene KXl lA. Schoter and H. J. Schulier 
;V Cu/r, Geoet 23. 375 (1993)). and the 'aerobic* 
, BOfOTTi of acetyl-CoA synthase. ACSl \M. A. van 
den Berg er a/.. J. Biol. Chem. 271 , 28953 (1 996)1 
with concomitant down-regulation of the glycoM- 

^ISf • ^ " • 5330 (1991)). Other genes 
7^^^^^^^ tn carbon metabolism but 
WKwn to be rtouceo upon nutrient limitation in- 
^nes encoding cyiosolic cataiase T CTTi 
^^ Bissmger et aJ..^ibiO. 9. 1309 (1989)1 and 
several gms encoding small heat-shock proteins 
' and HSP42 fl. Fan^V;; 

• SLioH^f^- 15602 (1991): U. M. 

S^^'^n^^ • ^^^^ Gen. Genef. 223. 

/ 2^^yK;^^^'^■•^•^•^-^^^ 

,37. The levels of induction we HDeasured for genes that 

^ icw levels In the uninouced 
^ f^®*r?^*^-^'«^'^')were9eneralfy|^ 
tt^^thwe presnously reported. This discrepancy 

w«ikely <kje to the consen^ative background sub- 
^ used, whicti generaity resulted r. 
ovwestmaioi of very low expression*^ 
38. Cross-hybndization of highty related sequences can 

occasionaly obscwe Changes in gene expres- 
aooanirTvortart concern vvhere of gene 

arefiJTK^^ specialized and drfferentiatty 

TS^^iJ^r^J!^ ctehyorogenase genes. 
^1 ana ADH2. share 88% nudeotide identity 

^^JH?^^ Of these genes is an important' 
feature of the diauxic shift, but was not obsen^in 
n« expenment. presumabty because ot cross-hy- 
onaratcn of the fluorescent cDNAs representing 
?i;?i!)^9«^-Ne^'ertt^ 
ttOdmerentiai expression of dosefy related isof orms 
ourtfier wzymes, such as hxki/hxk2 (77% iden- 

OAL7(73% identicaf) m. and PGM1/PCM2 (72% 
«*^;to|0|D Oh. J. E. Hopper. Mol. Ceff. So/. 10. 
1<l S (199Q)],n accord with previous studies. Use r 
deSberaiety selected DMA se- 
ouences corresponding to the most divergent seg- 

of homologous genes, in feu of the complete 
Oene sequences. Should relieve this problem in many 



40. a Tzamarias and K. StnrfH. Nature 369. 758 (1994) 

41. Orfterences in mRNA leveis between the tupiA and 
™:*yp^ strain were n^easured in two independent 
^^f^^ljents. The correlation coeffoe^ 
compiete sets of expression ratios measured »i 

these diofccate experiments was 0.83. The concor- 
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Miqf 

' hifl^ batwMn two 

o pan bam the atop oadanol««7lPi «M. ■ 

^^' ^a^^uJT^^' HV>tewa, MoL Gmi Genet 

.^.242.250(1994). , ' 

'^.?^S!5l!r A. 1^ M. j: Marines. ' 

V,n»l^'^ >W fi**oa Morofttt 60. 1783 . 

^I' ^ * 1 ». 3W (199lV 

4a jjLAWemmie. M. S. Szczypki 0. J. Thiete; W. s; 
J^je-Rowley. J. Bior Cham 269. 32592 (1094)r - 
^^^'^"^^ we scanned using' a^custom-tuit 

»aser.rniaoscope:bult by. S. Sn«h with ; 
.yft^;«gwnttan by N. gv. Dettfs co rriiiii> » «^ 
nod^andcoi^^ cmgm; 
. «antord.e<iVpbrowa acarvwd S'a 

-rwoMcn of 20 ♦im par pbiaL* A separate son. usm 

«wapp?Dpriaie excitation Ine, was dcna tor aaeft 5 
the two fluor o pliu a s r — • - 
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W the ratio between the signals h the two tf«v 

' !!^2u??L0«^ DNA. To normaftza the two 
^wnn«^ w^ respect to overaD Hen^ 
awed photomudipier and laser power aettras 

that the signal ratio at these tfOTwa wasS 
<*»e to 1 .0 as posstote. The oombrwd hwges wwa 
ana^y^J^ vvith ajstonr>-wrttten aoftvvM. A boi^^ 
box, fitted to the size of the Df^ spots h each 

ouadrartt^ was plaoed over aach »Tay tfemam. The 
average fluorescent intensity was ctfcUated by auFh 

fTwgma intensities 0* each pixel present in a bound- * 
«[»9box andthen dividing by the total nunter of 
poceis. Local area background was catoutatad 'tor 

each array element by determining the avwage auo- 
' 20% of ^ 

tow. AWxHjgh the method tends to mderestvnate 
tfw background, causing »i unoerestimatic^orS! 

treme ratios, it produces a very consBtmt and n«se- 
totewTt approximBtion. Although the analoo-to- 
d^rtal board used for data coiection poss^ a 
w»bB dynamic range (12 bUs). several signals were 

^^^.^^ Thera(ora. «eireiiS 

ratiM at bnght etemems are generaly underastrnat- 

It^^ ^ ^ «»«^^ in the 

oacKground measurements tor al elen«nts on ttw 

array. 

50. in ed^ to the 17 genes Shown ri Table 1 three 
additional genes were induced by an aver^ of 

rnore than threefold in the OuoUcate experimerttTbUt 

than twofold (range 1.6- to l.g-foW) 
"^^^^^V^ P- Speflman. J. Ravetto, M. 
^Ai^^®- ^ T, Feraa. and other mem. 
bars of the Brown tab lor their assistance and hetafii 
ad«a. We also thank S. Friend. D. Boiste^ 

smith. J.Hudson, and D.Dolginow lor adwioe auo- 
pwt^and ericouragemem: K. Stojhl md S. Owtti- 
lee tor the Tioi oeietion strain; L. Femanoes tor 
helpful advce on Yapl; and S. Ktaphotz and the 
rewewws for many help W comrnent^ 
scnpt. Supported by a grant (Torn the National Hu- 
Ger^ Research Institute (NHGRf) 

strtute (HHMI). J.O.R. was supponed by the HHM 
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insmuuonal Training Grant in Genome Science7Ta2 
HG00O44) from the NHGRI. pSHI^ISSS 
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Drug target validation and identification of secondary, drug 
" target effects using DNA friicroarrays " ' 
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We desi:ribe here a method for drug target validation and identirication of secondary drug tarV 
get effecu based on genome^de gene expression patterns. The method is demonstrated by 
several experimenu, induding treatment of yeast mutant strains defective in cal^ineurin inv 
munophilins or other genn with the immunosuppressants cyclosporin A or FKSOe Presence or 
absence of the characteristic drug 'signature' pattern of altered gene expression in drug-treated 
cells with a mutation in the gene encoding a putative target established whether that target was 
required to generate the drug signature. Drug dependent effects were seen In tarqetless' cells 
Showing that FK506 affeeu additional pathways independent of calcineurin and the im." 
munophilins. The described method permits the direct confirmation of drug targeu and recoo- 
nition of drug-dependent changes in gene expression that are modulated through pathwan 
distinct rom the drug's intended Urget. Such a method may prove useful in improvinrihe effi^ 
ciency of drug development programs. ' , ^i""?""- 



Good drugs are potent and specific: that is. they must have 
strong effects on a specific blolojgical pathway and minimal ef- 
fects on all other pathways. ponflrmaUon that a compound in- 
hibits the intended target (dnjg target validation) and the 
identification of undesirable secondary effects are among the 
main challenges in developing new drugs. Comprehensive 
nnethods that enable researchers to determine which genes or 
acUvities are afTected by a given drug might improve the effi- 
ciency of the drug discovery process by quickly identifying po- 
tential protein targets, or by accelerating the identification of 
compounds likely to be toxic. DNA microarray technology, 
which permits simultaneous measurement pf the expression 
levels of thousands of genes, provides a comprehensive frame- 
work to determine how a compound affects cellular metabolism 
and regulation on a . genomic scale' ". DNA microarrays that 
contain e»entially eveiy open reading frame (ORF) in the 
Saccharomyces cmvisiae genome have already been used success- 
fully to explore the changes in gene expression that accompany 
large changes in cellular metabolism or cell cycle progression 

In the modem drug discovery paradigm, which typically be- 
gins with the selection of a single molecular target, the ideal in- 
hibitory drug is one that inhibiu a single gene product so 
I completely and so specifically that it is as If the gene product 
were absent. Treating ceils with such a drug should induce 
changes in gene expression very similar to those resulting from 
deleting the gene encoding the drug s target. Here we have com- 
pared the genome-wide effects on gene expression that result 
from deletions of various genes in the budding yeast 5. cerevisiae 
to the effects on gene expression that result from treatment 
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With known inhibitors of those gene, products. Using the cal- 
cineurin signaling pathway as a model system, we tested an ap- 
proach thai permits idenUfication of genes that encode proteins 
specifically involved in pathways affected by a drug. The FK5d6 
characteristic pattern, or :signature*. of altered gene expression 
was not observed in mutant cells lacking, proteins inhibited by 
FK506 (for example, a calcineurin or FK506-binding-proteln 
mutant strain), but was observed in mutants deleted for ^genes 
in pathways unrelated to FK506 action (for example, a cy- 
clophllin mutant strain). Conversely, the cyclosporin A (CsA) 
signature was not observed in CsA-treated calcineurin or cy- 
clophilin mutant strains, but was seen in an FK506-binding.pro- 
leln mutant strain treated with CsA. The method also 
demonstrates that FK506. a clinically used immunosuppressant, 
has off-targef effects that are independent of its binding to im- 
munophilins. Thus, the approach we describe may provide a 
way to identify the pathways altered by a drug and to detect 
drug effects mediated through unintended targets. 

Null mutants phenocopy drug-treated cells on a genomic scale 
To test whether a null mutation in a drug target serves as a 
model of an ideal inhibitory drug, we examined the effects on 
gene expression associated with pharmacological or genetic in- 
hibition of calcineurin function. Calcineurin is a highly con- 
served calcium- and caimodulin-activated serine/threonine 
protein phosphatase implicated in diverse processes dependent 
on calcium signaling" In budding yeast, calcineurin is re- 
quired for intracellular ion homeostasis'*, for adaptation to pro- 
longed mating pheromone treatment'* and in the regulation of 
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5te«*unlt(erteine«lnA.«x«^ 

■Id calounvliinding regutatny submits cabnoduiin (CMD) and oidiwain B 
^)jM«r ««lng FK506;«k1 spedtaBy bind and.mhiblt the 
F?yy^*°^ tomerase actMty of thdr.napeeOve immunoptiBins. FKSOe' 

5*WorL Onig-inununophnn complens bind and hhM the calciunv and 
•"hr-djAv^Umulatrt^^ 

?neurtn are tramertpttottl actW^ that aa to niodulate gene expression. 



the onset of mltosl$'»./Jn mammals, caldneurin has been impli- 
cated In T-cell actJvaUon" In apoptosls". Jn cardiac hypertro- 
phy ...and m the transition from, short-tem to long-term 
memory". In both organisms', calcineurln activity is inhibited 
by FK506 and CsA. immunosuppressant drugs whose effecu on 
calcineurln are mediated through families of intracellular receo- 
tor proteiris caUed inununophillns'"' (Fig. j). To assess the 
fects of phaimacologlc inhibition of calcineurln. wild-type 5 
cwKisiae was grown to early logarithmic phase in the presence 
or absence of FK506 or CsA. bogenlc cells, from which the 
genes encoding the catalytic subunitt of caldneurin (CAM j and 
CAM2) had been deleted" (referred to as the cna or caicineurin 
mutant), were growri in parallel., m the absence of the dnia 
l^luorescently-labeled cDNA was prepared by reverse transcrij' 
Oon of polyA- RNA in the presence of Cy3- or CyS-deoxynu- 
cleotide triphosphates and then hybridized to a microarrav 
containing more than 6,000 DNA probes represenUng 97% of 
the Icnown or predicted ORFs in , the yeast genome 
Simultaneous hybridization of CyS-Iabeled cDNA from mock-" 
treated cells and Cy3-labeled cDNA from cells treated with 1 
Hg/nJ FK506 allowed the effect of drug treatment on mRNA lev- 
eU of each pRF .to be determined (Fig! 2a and A and data not 
shown) Similarly, elfecu of the caicineurin mutations on the 
mWMA levels of each gene were assessed by simultaneous hy- 
^"If !^'**'^k''' CyS-labeled cDNA from wild-type cells and Cy3- 
labeled cDNA from the caicineurin muunt strain(Fig 2d For 
„ each comparison of this kind, reponed expression ratios are the 
average of at least two hybridizations in which the Cy3 and Cv5 
fluors were reversed to remove biases that may be introduced bv 
gene-spedflc differences in incorporaUon of the two IIuor 
(dau not shown). 

Treatment with FKS06 in these growth conditions resulted in 
a signature pattern of altered gene expression in which mRNA 
evels of 36 ORFs changed by more than twofold 
(http://www.rosetta.org). A very similar panem of altered eene 
expression was observed when the caldneurin mutant suain 
WM compared to wlld-type cells. Comparison of the dianges in 
mRNA expression of each gene resulting from treatment of 
wl d-type cells with FK506 with mRNA expression changes re- 
suiting from deletion of the caicineurin genes showed the con- 
siderable similarity of the global transcript alterations in 
response to the two perturt>ations (Fig. 2M). QuantiHcaUon of 
this similarity using the correlation coefndent (p) showed 
large correlations between the FK506 treatment signature and 
the caldneurin deletion signature (p . 0.75 ± 0.03) as well as 
the CsA treatment signature (p - 0.94*0.02). but not wiUi a 
randomly seleaed deletion mutant strain (deleted for the 
YER07JC g^ne: p . -0.07 ± 0.04: Fig. 2e). The FK506 treatment 
signature was also compared with tiiose of more than 40 other 
deletion mutant strains or drug-treatments thought to affect 
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relations ?^ H^- ha" statistically s,gnincant<cor- 

cineunn function provides axlose and spedflc phenbcooV of 
treatment with FK506 or CsA. ^ pnmocopy or 

Jo^avoid generalizing from a single iiample. we also com- 

a«)le^(3.AT) .with U»e effects of deletion of the «J3 gene HISS 
en»des imidazoleglycerol phosphate dehydn-t^^S^ 

^^^IT. H "synthetic pathway In 

y^, 3-AT is a competitive inhibitor of tiUs enzymVti«t tirtg- 
^ ««criptional ammo-add starvatiwTtenjorise- 
Mia^array analysis of wild-type^. and isogenic Ws^eSSJm 
earns demonstrated the expeaed large genome-wlde tranSS- 
tion^l responses involving more than 1.000 ORFs) resS^ 
fi^m treatment witii a-AT (Fig. 3a) or from HISS deletion fflg 
afdt^niT'' """P^" °^^« 3-AT treatment slgnaSre 

0 76 . J?2")T«"'~''"" * •^8*' of corilation 

sm,ii *° 8enes that experienced 

mall changes in expression level (Fig. 3«. A, a negative cbritrel 
correlations between ti,e 3-AT treatment slgLture o^^e 
hb3 mutant signature and the caldneurin mut^t strain we^ 
not statistically significant (p - 0.09 * 0.06 and -O.ilTo 04^ 
spectiveJy). That both the caldneurin/FK506 and the hls3/i^ 
comparisons were highly correlated indicates tiiat in many «s« 

1 'H^"""'"" of Wlld-type celb treated with 7n 

inhibitor of that gene s produn. 



B^^^ "'««9y- Drug target validation with deletion mutanu 
Because pharmacological inhibition of different urgets might 
g ye Similar or identical expression profiles, simple cSr^pS 
of dnig signatures to mutant signatures is unlikely to unambiTu 

addtina?? "^""'"^ this' imitation 'an 

additiona decoder step is used. We fim compare the expres- 

s on profile of wlld-type drug-treated cells to the expresltoj p^ 

nies from a panel of genetic mutant strains, using a conelai^S^ 

coefncient metric. Mutant strains whose expression p"onS 

sZl ^ . H °' "'"^■"'''''^ are selected anS 

subjected to drug treatment, generating the dnig signature in 
the mutant strain (that is. the mutant drug signature) If th^ 
muuted gene encodes a protein Involved In a pa'thway aff^^S 
by the drug, we expect the drug signature In mutant celU to be 

sin?"' M ^""^ ««~« Signature 

seen in wild-type cells. 'B'wiure 
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fiQ* 2 • Eipression profiles from 

fK506-trmed wild-type - (wt) 

celb and a catctneurin-disruptJon 

muunt strain ihare a genoim* 

widecorreiatioa' DNA miaoerray 

•analysis showing changes in gene 

expression resulUng from FK506 

treatment-(« and or from .ge- 

netic disruption of genes encod- 
ing . cateineurin «. - Pseudo- 
^color image of the results of si*. 

muttaneoiB hytiridteation of Cy5- 
'tat>eled 'cDNA (red) from . 

mocli-treated strain. • R563 and , Cy3-lat>eled xONA 
: (green) from strain RS6d treated with.1 . iig/ml FK506. 
.ib, Enlarged view of the boxed area in e. Arrowheads in- 
dicate specific ORFs induced or repressed, e. Pseudo- 
color image of the results of simultaneous hybridization 
of CyS*labeled cDNA (red) from s^^^in R563 and Cy3- 
tabeled cDNA (green) from strain MCY300 (deleted for 
the CNA1XNA2 catalytic subunfts ' of calcineurin}. 
Arrows indicate specific ORFs induced or repressed, 
'The iogw of the expression ratio for each ORF derived 
from the FK506 treatment hyt>ridtzatjons is plotted ver- 
sus the logw of the expression ratio in the calcirfeurin 
mutant hybridizations. ORFs that were induced or re*, 
pressed in both experimenu are shown as green arKl. 
red dots, respectively, c The log„ of the expression ratio for each ORF de- 
rived from the FKS06 treatment hybridizatioru is plotted versus the log„ 



wt VL calcinmtninuM 




E 



t 



r 




2 



0 - i 

1 ■!- 




l.og« dUC) calcineurin muution 



tog* (R/G) Jcmuution 



Of the expression ratio in the >ef077c rTiuiim:hybridi2aUons No 
were induced or repressed in both experiments. 



ORFs 



To Illustrate this, we. treated the Jjis3 mutant strain with 3- 
AT. The signature pattern of altered gene expression resulting 
from treatment of the mutant strain with 3-AT was much less 
complex than that of the 3-AT signature in wild-type cells (fig. 
4). This is seen simply by examining plots of mean intensity of 
the hybridiratlon signal (which approximately rcflecu level of 
expression) versus the expression ratio for each ORF (Fig.. 4). 
Genes that were expressed at higher or lower levels in 3-AT 
treated cells or in hJs3 mutant cells are shown as red and green 
dots, respectively. We analyzed the 3-AT signature in wild-type 
(Fig. 4a) and hJs3 mutant cells (Fig. 4c). as well as the bis3 mu- 
tant strain signature (Fig. 4b}. Whereas histldine llmitaUon in- 
duced by 3-AT induced more than 1,000 transcription-leve] 
changes in the wild-type strain, few or no transcript level 
changes were induced by treatment of the /iisJ-deleUon strain 
with 3-AT. This indicates that with the growth condiUons used, 
essentially all of the effects of 3-AT depend on or are mediated 
through the HIS3 gene product. 

Applying this approach to the calcineurin signaling pathway 
showed the specificity of the method. The calcineurin mutant 
strain and strains with deJeiions in the genes encoding the 
most abundant immunophilins in yeast*' {CPHJ and FPRJ) 
were treated with either FK506 or CsA to determine the profiles 

Table 1 Signature correlation of expression ratios as a result of FK506 

treatment in various mutant strains 



wild-type 
<ft/-FK506 



ene 
4/-FKS06 



wild-type 
FK506 



fpr7 
4/-FK506 



enafpri 
♦/-FK506 



0.93 1 0.04 ^.01 1 0.07 -0.23 i 0.07 0.12 i 0.07 0.79 i 0.03 



*^ *r!l^ °' "8^'"^' »pecif«lly m ,ne ailc.neur,n Una) and tprl 

(major FK506 bmdtfig protein) MMn muttnii. cnan^nu the mutem wim deletions of the caiaM-c wb 
unju Of ttjciheurin. CNA1^CNA2, The correuton cocffident repoaed in the first column re^e«n«?r^ ^c^ 
reUtKW between tvK) pain of hytyidi2at*om from inOepenoem wild-iype W- FKS06 eiperimenu 
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of altered gene expression resulting from drug treatment of the 
mutant cells (that is. mutant W- drug). We compared the drug 
signatures in the mutants to the wild-type drug signature using 
the correlation coefficient metric (Table 1). Although the signa- 
ture generated by treatment of wild-type cells with FK506 was 
highly correlated to the calcineurin mutant strain signature (p 
- 0.75 ± 0.03). it bore no similarity to the profile after treat- 
ment of the calcineurin mutant strain with FK506 (p » -0 01 i 
0.07). This indicates that FK506 was unable to elicit its normal 
transcriptional response in the calcineurin mutant strain 
Likewise, treatment of the fpr} mutant strain with FK506 
elicited an expression profile that was not con-elated to the 
FK506 signature in the wild-type strain (p - -0.23 ± 0.07), indi- 
cating that the FPR I gene product is likely to be involved in the 
pathway affected by FK506. The same was true for the cna ij^rJ 
mutant strain. In contrast, treatment of the cphJ mutant strain 
with FK506 generated an expression profile highly correlated 
with the wild-type FK506 expression profile (p « 0.79 ± 0.03), 
indicating the cpbJ mutation did not block the mode of acUon 
of FK506 and thus is not directly involved in the pathway af- 
fected by FK506. We tabulated the change in expression in re- 
sponse to FK506 in different mutant strains for all ORFs with 
expression ratios greater than 1.8 in FKSOe-ireated cells or in 

the calcineurin mutant strain (Fig. 5a).The 
calcineurin mutant strain signature and the 
FK506 responses in wild-type and the cphl 
mutant strain are similar, and there are no 
transcript-level changes (seen in black) for 
treatment of the calcineurin. fpri and cna 
fprl mutant strains with FK506 (Fig. 5a). 

Similar experiments and analyses with CsA 
provided further validation of this approach. 
The expression profile elicited by treatment 
of wild-type cells with CsA was highly corre- 



cphl 
-^/-FK506 
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Fig/ S 'Expresskm phifHes 
froiri -J /liEsi ' mutant iuiln 
•nd ^ wiw-iype . oils 
tretttd wmi share : « 
genome-wide rcorretaticm. 
OnA : microiirray ' analysis 
showing ichanglurm i gene 
exprwion resulting from 3*; 
AT treatment (a) or from ge- 
netic disruption of the Hi$3 
gena . (c). e« Pseudocolor 
image of the resutts of simul- 
taneous ' hybridization of 

CyS-labeled cDNA (red) from mocii-treaied wild-t^ suain R49i and 
Cy3.labeled cDNA (green) from strain R491 treated with 10 mM 3.AT 
*, Plot of the log« of the expression raUo for each ORF derived from the 
3.AT ueatment hybridizations is plotted versus the log« of the expression 
raUo in the his3 mutant hybridizations. ORf s that were induced or re- 
pressed in both experiments are shown as green and red dou resocc 
tively. The correlaUon of expression ratios applies not only to genes wrth 
large expression ratios (for example. CAM 7 and ARC 7). but atso extends to 
genes with expression ratios less than 2 (for example, tm and CPH1\ 
tm is induced 1.9.fold and 1 .5-fold, and CPHl h downregulated 1 9.fo!d 




-■J J.' -J. 



•4 . 



WIVL 




tna 1 7.foW, in celb ueated with 3.AT and AU i^m «us. resoeciiviriv 
T«« 0RF5 do n« fan on tt« Mne x - y. The l.ftmo« polmTJ^SSS 
P«m. Which i, induced by 3.AT ueatmem b« J^bn^^^ 

•» ^°«03w. Both dtt.^« 

m^Nl f^n^r"^"'"""""' ^0«03.«»mostii.*y^«^ 
mRNA. as YOR203W overlaps the H/J3 open reading tntZVp^Ju^ 

color image of the resulu of sirnulUneomhybridiXJ^bilTtS^ 
DNA (red) from wild-type strain R491 .nd1^3.2S^?D?MaSS 
from strain R1226. deleted for the H,53 gene. Lrowheads SaK 
ciHc ORFs induced or repressed. "wncaie spe< 



ated to the pronie elicited by miitaUon of the caldneurin genes 
0> - 0.71 ± 0.04). but did not correlate with the expression pro- 
flie resulting from treatment of the caldneurin mutant strain 
With CsA (p . -0.05 . 0.07: Table 2). indicating that the ^e^tl; 
deletion of caldneurin interfered with the ability of CsA to 
elidt its normal transcriptional response. Likewise, the CsA sic 
nature was essentially absent in CsA-treated cphl mutant cells 
and the expression proHle of CsA-treated <phj muunt cells cor- 
related poorly to that of CsA-treated wild-type cells (p - 0 18 ± 
0.07). Thus, the CPHl gene product was required for the CsA re- 
sponse seen in wild-type cells. Conversely, treatment of fprl 
mutant cells with CsA resulted in an expression pattern verv 
similar to the profile of CsA-treated wild-type cells (p - o 77 i 
0.03). indicaUng that FPU! was not necessary for the CsA-medi- 
aied effects. Analysis of individual ORFs affected by CsA and 
their expression ratios over the entire set of experiments con- 
firmed that CPHl and the genes encoding caldneurin. but not 



FPR]. are nec^o' for the wlld-type CsA response (Fig. Sfc). The 
observation that the profiles resulUng from FK506 or CsA drug 
reatment are similar to that of the caldneudn deletion mutant 
tram might allow the predicUon that caldneurin was invol^ 
n the pathway affected by these drugs. But because the expres- 
sion profile of the fprj mutant strain did not bear a stmng siml- 
lariy to the wild-type dnig expression profile for FK506. It is 
obvious that the drug treatment of the mutant strains was nec- 

~ \° 'l!""'^ '^P*''- " FK506 drug 

targe m the same way. the decoder" strategy was neccssaiy to 
.denUfy Cphl. but not Fprl. as a potenUal dTug urget for CsA. 

•Decoder- approach can identify secondary drug effecu 
For a drug that has a single biochemical target, the strategy out- 
lined above may be useful in target validation. In manyaaes 
however, a compound may affect multiple pathways and elidt 
a veo, complex signature. Decoding' such a complex signanw 



wt •/« lOmM 3-AT 



:£ 
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.2 -» 



wi vs. fto3 mutant 



his3 muunt 10 mM 3^T 




Log. (intensity) 




nIL . .??r* ^ """""^ nearly com- 

tZI^ V * P»°^ «^ ^ '<>9h Of the mean intensity of hy 

'"S'^ ^''P'^"*^" ^or each 
rl^^ J" "^"^ *° * pseudocolor image of a representative 
portion of the microarray. ORFs that are induced or repressed at the 95% 
conhdence level are shown in green and red. respectively, a. Expression 
''n.?!^ f;«""/'e«^nt of the wild-type (wt) strain with 3.AT. CyS labeled 
cDNA (red) from mock-ueated strain R491 and CyS-labeted cDNA 
(green) from suain R491 treated with 10 mM 3.AT. t. Expression profile 
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Logw (intensity) 



Log. (intensity) 



from the /./S3 deletion strain. CyS-labeled eONA (red) from suain R491 
and Cy3.lat>e.ed cDNA (green) from su.in R1226. deleted tTZ nil 
gene. . Expression profile of treatment of the hsS deletion strain with 3- 

1 "''^l^e'e- "rain R1226 nid Ss J- 

beled CDNA (green) from strain R1226 treated with 10 mM 3.AT 

TK ""^ corresponding to 

rors tend to inoease rapidly because spot intensities are not sufficiently 
above t>ackground intensity. »«"iicierniy 
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Table 2 Signature correlation of expression ratios as a result of cIJ 
treatment in variotB mutant strains 



wUcMype 

*/-ClA 



«/-CsA 



oia 
♦/-C»A 



tptl 
«/-CsA 



cnatphi 
♦/-CsA 



0.94 « aCM ^.OS 1 .07 0.77 1 0.03 



ephi 



-o n « 0.07 0.18 1 0.07 



''9'Wiirt camlitianihoimuieabicnt* Of the CiA <■«... ' 

«»«»««n ««.o pUnof tyt>iflmtlon» from tnotptnOg^^^'^J^J^^!^'"^ "* cofreliiion 



Into the effects mediated through the Intended tarcet (the on 
urget Signature-) and those mediated through unimended t^r 
gets (the off-targef signature) might be useful in evarualra 
compound's specificity. Our decoder' strategy is based on thp 
premise that off-targef signature should befnserTtive « tSe 
genetic disruption of the primary target. 

To detemune whether the decoder- approach could idcntlfv 
an oir-targef profile, we looked for a'^Lg-responsive S 
whose expression is Insensitive to deletion of the pZanf tar 
get. To increase the likeUhood of observing such ge" e7 he 
same strain described in Tables 1 and 2 were treated w^J, 
Wgher concentrations (50 Mg/ml) of FK506. This led to a mTch 

tL t at'Z^^h"'""'''" "^'-'y^' ''"'^""ng 

!iinp iiS^^* r """T"""- ^ ^"Wbiting or acti 

S .nH H °^ 'his expanded 

FKS06-induced expression profile were not affected by the cal 
dneurln. cphj or fpr, mutations, as drug treatment ofThlse mu: 
unt strains did not block their presence in the fS 
expression signature (Fig. 6). This Indicates that FK506 was u,/ 
gering changes in transcript levels of many genes throueh oath 
ways independent of calcineurln. CPHl Bn! FPR, MTlft: 
up«gulated ORfs in ti,e off-targef pathway we« gene 'r' 

3 S^f 2,)'", ^ transcriptional activator Ccn4 

(ref. 24). In «,me strains, a reporter gene under CCN4 control 
was induced in response to FKS06 treatment" To de,.™- ° 
whether CCN4 Is involved in this pathway SLt Is^ndeSm 
of caldneurin. CPH, and FPRl. we analyz^ U,e SfecJof u«T 
ment with high-dose FK506 on globd^ »ne exon^i^n 
strain wlU, a CCN. deletion (Fig.?,. Of t?e 4 1 STv^^h cal' 
cineurin-independent expression ratios greater than 4 32 vSe 
not induced in the ^4 mutant, indicatinp that t hJ^r inH T 
by FK506 was CCAT^ependent. Not aTic^S^S "^^^^ 
were Induced by FK506. This FKSOe-induced ,ulJt of Ccn7 

FKSOfi l °^P"*«P» "her regulatory circuits preven 
FK506 activation of some CCN^-regulated aenes Sev,.n nf .i? 
remaining nine ORFs induced by «506 wf^TndtTendem oJ 
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both the caldneurin and CCN4 pathway. The 
simplest explanation is that FK506 inhibits or 
activates additional pathways. Members of thb 
class include SNQ2 and PDXS. genes that «^ 
code dnig efflux pumps with structural homoJ- 
ogy to mammalian multiple drug resistance 
proteins". FK506 may interatt directiy wlti, 
PdrS to inhibit its function". Our results indi- 
cate Uiat treatment with FK506 leads to four- 
fold-to-sixfold induction of PDWmRNA levels 
rORl. another gene tiiat can confer drug reals- 

newly discnLr«r >"c'"*n« the nature and extent of 

SdiS thH;;;."'' ""^""''^ P-'^'-y' ^■ 

iZlrZ?*^ ^'^ ' ""'S '"'Ket validation and the 

Wentification of secondary drug target effects that uses DN A nS 
CToarrays to survey the efTecu of drugs on gloW «ne 
sjon panerns. We established thar oot-.., "o" gene expres- 

inhibition Of gene runcS'o? c^ rLT tS;:eTyT:r 
changes in gene expression. We also demonstr««S?h« 

confirm a potential drug target by treTnTrrf!?. . 

dpfprtii/* i« rK* "wgei oy ireating a deletion mutant 

at^ Inli P"""^« "'set. Drug-mi" 

ated signatures from strains with mutations in pathwavs or 
processes directly or indirectly affected by the dn.gC mJe Z 



S 



•c 

E 
< 



to 
Z 

00 
Oi 
Oi 



Q Stntn: 
FKS06: 

'»JH7IC 

TP»» 
TPSi 
YBRmC 

nni»«c 

PftBi 
CA«C2 

vuAattw 

TBAOOSW 

L Strtrn: 



cna 



wt 



cphi 



cna cnafprl 




cna 



wt 



cphl 



cn» cnafpri 



green (indunion)-,ed (repression) t^S^^r^ ^ 

FK506 signature genes have expresion ratios near unitv inZZ^TT^ 
involved in p„hw.„ affected by FKSOe ic^J^. J^I^ 
tants) bu, in deletion su««1^ unrelau^S^S,ti,7.'^*r ' 
(cna) mutant and CsA treatment signau^rS^i^ta L r " 
columns. Almost .« CsA «g«,.ure gen^ r«ve eS^rraL^"ll"! 
delet«n susin, involved in pathway, affeaed t^sA (calc^r^,? - 
cna muu.™,, but no. in deletion suains ^Ji^^^S^^^ 
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no similarny to the wild-type drug expression profile -in con 
trast. drug-medlated signatures from strains with mulaUo,« in 
genes involved In pathways unrelated to the d^gf aSon 
showed extensive similarity to the wild-type drug sisnatu^ 

(FK506). we were able to decode a complex sigHature Into coiT 
ponent parts, including the IdentfflcaUon of an on-^ZS 
nature that was mediated through pathways independent of 
caldneurin or the Fprl Iminunophlliri.' ."^Penaent of 

Discussion 

It is well-established that high-throughput biochemical screen 

target. The decoder' approach described here complement 

spedflc^ty: the tendency of a compound to inhibit patlTwaJ 
other than that ofju intended target. The ability to oSe 
such off-target- effect, will likely be useful in several ways 
Pronimg compounds with known toxicities will allow the de 
velopment of a database of expression changes associated with 

ofT-target signatures of othenvise promising compound then 
may allow earlier Identification of those likefy to fall ?n djnlca" 
trials. Comparing the extent and peculiarttie/of off-tarRet 1 
natures of promising drug candiates could provide a nf w way 
to group compounds by their effects on secondao' pathwrys 
even before those effects are understood. This may p^veTo be 
an altemathre. potentially more effealve. way to select com 
pounds for animal and dinlcal trials. Some drugs are r^o^Tf ' 
ective against a related protein than against the orT^y 
intended target. Sildenafil (Viagra™), for example, was fSll 
developed as a phosphodiesterase inhibitor to control ciS« 
contractility, but was found to be highly specific for phospho 
diesterase S. an isozyme whose inhibition overcome, !efecu in 
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of 4 in at lean one e»p«ime« w Used and th*«n^^ 

indicted s,«i„ are *own in J^SSS!!!^ 
«on) cotor scale. The genes have beST^STSTSS^ 
spending to thew expected beha»ion.-^^Otf*I^fS? 

to FK506 (50 H9/mO except wh^iefSSSSToeTS 
/W» Of both •re tWeied; •CCW4KJependem^o«« m^.!^ 
ex^pM^^GOV. ^^eietedTSSL^rSTt??^ 
wt*n Mle.ne.*», genes orTW, or CPH1 m deleUrtihTt 



penile erection. It is possible that application of the de 
coder to other compounds may shJ. that th^y tooLe a 

t^entdtrrgT ^«^^"« .' .-«« 

enr f.'--^™^^^ 

.quence-www.,tt„ford.edu/group/yeMt.deleUon.pS 
ject/deieuo„.html). Efforts are underway ^obtain 
expression profile,, from each deletion mut^t It^n 

^n^o " '""'^"S Inactivaa^T^: 

senual genes presents a unique problem but It ^ 

possible to do so by examining hetlroz^gSSTr' by^^a'^oS^ 

Ihh' " «P«"'°" of «he e^entS gew 

AJthough „ „ already feasible to test several com^unSrin 
dozens of yeast strains, another challenge for the •dJ^ode," 
tra.egy.wlll be the effident selection of th? muuno wS^SS 
oons .n genes most likely to encode the Intended d-^g SiSt 
The signature con-elaUon plots described .» J/ T 
could be used aspart of thrt select';?5^ocS ZZl" "^! 

sents additional challenges. It is considerably more dim2uh?„ 
«plate functionally targeUess" cells. Strategi^rvolZ 
able promoters, known specific inhibitors. anti-semrSrs ri 
bozymes. and methods of tareetino .nJ^L . , " 
degradation are possible and shoSd tSS LoTr .Tm. " 
tjon » that not all cell types express the same stTgen^^nd 

tJ. fin h1 """y d'fferJm cell 

types. In addition, applying the decoder' to human cdh i» 

aUo require technical improvement, that allow expr«„on pT! 
f. hng from a small number of cell,. Even the broa?eTQuestr„ 
of whether the insensitivity of on^-targef slgnVtuJ« o .hTJ? 
ruption of the main target i, the excepUon or "Se c J^^^^^^^ 
be answered by the accumulation of more dau BarkT, ° 
Leibler. however, have argued in favor of roUness oVSi.::, 
cal network,, indicating that drug perturbationr( off « !f : 
signatures) may be robust even when the system is sub ecS^L 
another perturbation (such as a genetic dlsrup t oXf 28) 
Many practical developments will be neces«,y if the ScodI ^ 
concept I, to be broadly applied. oecoder 
Expression arrays have been used mainly as an initial screen 
or genes induced in a particular tissue or process o? in Ij^bv 
ocusing on gene, with large expression ratios We haS 

nd"reoearer' '° experimental p^toc^ 

and repeat experiments increases the reliability of the data and 
permit, new applications. For example. It proWdes a l^ger 
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R563 
R558 " 
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MCY300n • 

R132 ^ 
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BY4719 
6Y4738 
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BY4728 T :j : 
8X4729^ K 
R1226 i 



Mevant genotype 

^ A4ato.i^5.5P /)4?i507 «W-707 l^»7:itej AdJrilPOO fcu^^'f 
Atetf un3:52 1^2-801 mk2-1p1 tipiJ^ his3'A200 ieu2'^ 7 /ttij-H/SJ 
A4Ma ir^J.52(|«^^7;«d^.7p7^t^ 
; . Mats ^^'S2lp2;8madt2'10ltn>h^m^ . , 

Mata trpl'^ 63 ura3^ \ : ^ ; s ' s 

M0ta/aBY47l9JrBY4738 . 
Mata his3'A200tTpl'A63 ura3'A0 
Mata hi$3-^ 20Q,tfplTA63 un3'A0 
Mata/aBtA72B XBU729 . 
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of genes at higher confidence levels ' that serve as a more 
linlque signature for a given protein penurbaiion. In addiUon. 
it allows subtle signatures to be detected, when, for example, a 
protein is orily partially inhibited. This may enable clinidal 
monitoring of small changes in protein function in disease or 
toxicity states before they could otherwise be detected. 
Because the functions of many genes detected on transcript ar- 
rays are known, these mi croarrays are powerful tools that pro- 
vide detailed information about a cell s physiology. For 
example, changes in the flux through a metabolic pathway are 
renccied in transcriptional changes in genes in the pathway'. 
Furthermore, it may be possible to indirectly measure protein 
acUvlty levels from expression profiling data (S.F.. a/., un- 
published data). Thus, although the eventual development of 
genomic methods allowing the direa measurement of all eel- 
lular protein levels will be an important achievement, tran- 
script array technology offers an immediate and robust means 
of evaluaUhg the effects of various treatments on gene expres- 
sion and protein function. 



ATv « indicated. Celh were harvested and « mRNAJsolaied as at>ove. 
FKSOe was Obtained from, the Swedish ; Hospital Pharmacy = (Seattle. 
W«hingion) end purified io homogeneity by ethyt acetateeSraction^^ 
I Simon (Fred Hutchinson Cancer Research Center. Seattle. Washington) 
CsA.was Obtained Irom^Alexis Biochemicals (San Dieoo. California)- 3^T 
was from Sigma. • 'l*. 



1","!!"!^. "y^ri"""'*" o' the tainted Mrnpl.'. FluorescemiH.. 
betedxDNA was prepared, purified and hybridized esientially aj de- 
scribed'. Cy3. or CyS^UTP (Amersham) w« incorporated into cONA 
durmg revene transcription (Superscript Jl: Ufe Technoiofliei) and puri- 
f«d by concentrating to less than.lO mI using Mlaocon,30 microconcL 
Houston. Texas). Paired cONAs were resus;K«ded in 
20-26 Ml hybr,dttat.on solution (3 ?< SSC. 0.7S pg/ml poiyA DNA, o:2% 
h.f « 5 """"""y "n«»er a 22- x SO-mm coverslip ror 6 
nat63 Call according to a published method'. ■ 



K Methods 

Consmjction, growth and drug Ueatment of yeast strains. The strains 
used in this study aable 3) were constructed by sundard techniques" 
To consioia strain R559. strain R563 was transformed to Leu- with plas- 
mid pM12 digested by SaA and MM (provided by A. Hinnebusch and T 
Dever) . Strains R1 32 and Rl 33 were constniaed by Uansf or ming the bac- 
terial Itanamycin resistance cassette" flanlted by genomic DNA from the 
CPHl and FPRl lod. respectively, and selecting for G4l8-resisi8nt 
colonies. For experimenu with FK506. cells were grown for three genera- 
lions to a density of 1 x 10' cells/ml in YARD medium (YPD plus 0,004% 
adenine) supplemented with 10 mM calcium chloride as described'* 
Where indicated. FK506 was added to a final concentration of 1 jig/ml 
0.5 h after inoculation of the culture or to 50 pg/ml 1 h before cells were 
collected. CsA was used at a final concenuaiion of 50 ug/ml Cells were 
broken by standard procedures" with the following modifications: Cell 
pelleu were resuspended in breaking buffer (0.2 M Tris HCI pH 7.6. 0 5 M 
NaCI, 10 mM EOTA. 1% SOS), vortexed for 2 mtn on a VWR multi-iube 
vonexer at setUng 6 in the presence of 60% glass beads (425-600 pm 
mesh: Sigma) and phenol.chloroform (50:50. volume/volume). After sep- 
aration of the phases, the aqueous phase was re-extracted and ethanol- 
precipitated. Poly A* RNA was isolated by two sequential 
chromatographic purifications over oligo dT cellulose (New England 
Biolabs, Beverly. Massachusetts) using established protocols". 

For experiments using S-AT. wild-type or his3/his3 cells were grown to 
early logarithmic phase in SC medium, pelleted and resuspended In SC 
medium lacking histidine for 1 hr in the presence or absence of 10 mM 3- 
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corn^n 'nicroarr.ys; PCR products conuining 

common 5 and 3 sequences (Research Genetics. Hunisville. Alabam^ 
were used as templates with emino-modified ton<vard primer and unmod- 
ified reverse primers to PCR. amplify 6.065 ORFs from the S. cetB^siae 
genome. Our firsi-pass success rate was 94%. Amplification reactions that 
gave products of unexpected sizes were excluded from subsequent analv* 
SKS^ORFs that could not.be amplified from purchased templates were am- 
plified from genomic DNA. DNA samples from lOO-pl reactions were 
isopropanol-precipiiaied. resuspended in water/brought to a final con- 
centraiion of 3x SSC in a total volume of lS pt, and transferred to 384- 
weM m^oiiier plates (Genetix Limited. Christchurch. Dorset. England). 
PCR products were spotted onto 1 x B-inch polylysine-treated gtasi slides 
by a robot built essentially according to defined specificaUons**' 
(http://cm9m.stanford.edu/pbrown/MGuide). After being printed slides 
were processed according to published protocols'. 

Microarrays were imaged on a prototype multi-frame CCD camera in 
development at Applied Precision (Issaquah. Washington). Each CCD 
image frame was approximately 2.mm square. Exposure times of 2 s in 
the Cy5 channel (white light through Chroma 618-648 nm excitation fil. 

TrLnTT^^^!!:'^^'^ ""^ ^''^'^ '"^ 1 s in the Cy3 channel 

(Chroma 535-560 nm excitation filter. Chroma 570-620 nm emission fil- 
ler) were done consecutively in each frame before moving to the next 
spatially contiguous frame. Color isolation between the Cy3 and Cy5 
Channels was about 100:1 or better. Frames were 'knitted* together In 
software to make the complete images. The intensity of spots (about 100 
nm) were quantified from the lO-pm pixels by frame-by-frame back, 
ground subtraction and intensity averaging in each channel. Dynamic 
range of the resulting spot intensities was typically a ratio of 1 000 be- 
iween the brightest spots and the background-subuacied additive error 
level. Normalization between the channels was accomplished by normal- 
izing each channel to the mean intensities of all genes. This procedure Is 
nearly equivalent to normalization between channels using the intensity 
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filio Of genomic DNA spots', M h possibly more robust, as it is based on 
the biterttUies of several tfrausei^spoudQtributedom them 



Stgnature. correlation cbefrtcienti^.and thifr^eor^flctew linto. 
Cdrrelation'coerncieftts between the; signature ORFs of v^ious experi^ 
ntertts were catculatod using:, f-^. . ' - '.^i,^- ' 

where is the log« of the expressionVatio for the k* gerw in the x signa- 
tuo, and is the tog^ of the expressidn ratio for the k-tgene in the y sig. 
nature. The sumrriation is dyer those genes that were either up- 6r' 
down-regulated in either experiment at the 95% confidence level These 
genes eachhad:a less than 5% chance of being andatly unregulated (hav- 
ing expression ratios departing from unity diie to measurement errors 
•tone). This confidence tevet wAs assigned based on an error model which 
assigns e tbgnormal probability dttiribuiion td each gene's expression 
ratio w*th,char8Cteristic width based on the bbserved scatter in its re- 
peated measuremenu (repeated arrays at the same>iominal expeftmental 
conditions) and on th^ individual array hybridizetiori quality. This latter 
dependence was derived from control experiments in which both Cy3 
and Cy5 samples were derived from the same RNA sample. For large 
numbers of repeated measurements the error reduces to the observed 
scauer. For a single measurement the error is based^on the array quaiitv 
and the spot intensity. • . . ... 
. Random measurement errors in the x and y signatures tend to bias the 
correlation towards 2ero. In most experiments, most genes are not signif- 
icanUy affeaed but do show small random measurement errors. Seleaing 
only the '95% confidence' genes for the correlation calculation rather 
than the entire genome, reduces this bias and makes the actual biological 
correlations more apparent. 

Correlaiiom between a profile and itself are unity by definition. Error 
limits on the conelaiion are 95% confidence lirniis based on the Individ- 
ual measurement error bars, and assuming unconelated errors" They do 
not include the bias mentioned above: thus, a depaaure of p from uriity 
does not necessarily mean that the underlying biological correlation is ihv 
perfea. However, a correlation of 0.7 ± 0.1, for example, is very signifi- 
canlly different from zero. Small (magnitude of.p < 0.2) but formally 
significant correlation in the tables and text probably are due to small sys- 
tematic biases in the Cy5/Cy3 ratios that violate the.assumption of Inde. 
pendent measurement errors used to generate the 95% confidence 
limits. Therefore, these small correlaUon values should be treated as not 
significant. A likely source of uncorrected systematic bias is the partially 
correaed scanner detector nonlinearrty that differently affects the Cv3 
and CyS detection channels. . 

The 1 |ig/ml FK506 Ueatment signature was compared with more 
than 40 unrelated deletion mutant strain or drug signatures. These con- 
trol profiles had correlation coefficients with the FKS06 profile that were 
distributed around, zero (mean p . -0,03) with a standard deviation of 
0.16 (data not shown), and none had correlations. greater than p . 0 38 
Similarly, the calcineurin mutant suain signature correlated well with the 
CsA treatmem signature (p . 0.71 ± 0.04) but not with the signatures 
from the negative controls (mean p . -0.02 with a standard deviation of 
0.18). 



wtwlter spots have fewer Image pixels in the .vefiBe. ThH do« not ite. 
grade eeeuracy noUca6ly untR the number oT Wb bitb* 
wh«:h«,e the spot tnwn the diu leTTW-Sj^SIw 

^l^^r^^ " the nomine grt«Jii ,d.pUv^ kicked ^ ^^X' 
regions by the image processing »ttwai.r jUriequrt .pot -weiider- «Ah^ 
l^bregion greater than h.«:..,pot seeing b.a dlBleuny ter t^Jw 

•na^.based on human Inspeaion oTfte 'winder-. Any sJoTpertwI 
ov?rtopp«^ are excluded from.the data set Less than ! »x7»potti^ 
ca Ware »ipptiui fnr Ik... ' . ■. ■> .• •t'^.^jp*- 



cally'are rejected, (or these reasonsi 
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Quality conuots. End-to^nd checks on expression raUo measur^nt 
accuracy were provided by analyzing the variance in repeated hybridiza- 
tions using the same mRNA labeled with both Cy3 and CyS and also 
using Cy3 and CyS mRNA samples isolated from independent cultures of 
the same nominal suain and conditions. Biases undetected with this pro- 
cedure. such as gene-specific biases presumably due to differential incor- 
poraiion of Cy3. and CyS-dUTP into cONA. were minimized by doing 
hybridizauons in fiuor-reversed pairs, in which the Cy3/Cy5 labeling of 
the biological conditions was reversed in one experiment with respect to 
the other. The expression ratio for each gene is then the ratio of ratios be- 
iween the two experiments in the pair. Other biases are removed by algo. 
rithmic numerical de-trending. The magnitude of these biases in the 
absence of de-trending and fiuor reversal is typically about 30% in the 
ratio, but may be as high as twofold for some ORFs. 
Expression ratios are based on mean intensities over each spot. Some 
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CO moMic WW RNA WIS obuintd ly phenol and 
cMorofonn manaion of the vims and preo^tcd 
fromcthMoLCA^uscmbty rasaiDni in tht pro- 
•net of nonccffnitt RNAs were identical to those 
f »^ in (9). in the abscnct of WW. CA^ conts 
formed under the foUowinf contftiom: 300 »iM CA- 
NC 1 M Naa and SOmM tri»-«a (pH ao) at 37X 
tef 60 min. in the absence of nogcnoui RNA, neither 
cones nor cylinden formed at concentratiora of OS 
M N»a or below. Absorption spectra demonstrated 
that our CA^C preparations were not contaminated 
with CscMcMt coff RNA (estimated lower detection 
tonh was *1 base/protein molecule). To control for 
•ven tower leveU of RNA contamination, we prein- 
cubated the CA^ protein with 0 J n«M4 ribonu- 
dease A (Type VAS, S4 Kimitz U/n^ Sl«ma) for 1 
hour at 4«c which then formed cones normally 
13. V. Y. icttshko. dau not shown. 

M. Ce and K. Sattier. Om. Phys, Utt, 220 192 
(1994). 

15. A. Krishnan er »L. N»tun 383. 451 (1997). 

16. L B. Kent ft aZ../ Mrot 7Z 4403 (1998). 



n*^ Wotted bnefiy with filter p,p„. ^ 
bqu.d jrthane. and transfmed to liquid liZk ™ 
xen fndi were transferred to a Philips 420 TtM 
•WPJ^ « Catan coldiuie syitenx and images 
of p»ndtt « wreous ice were recorded under tow 
dose conditions at 36.000X magnification and -1^6- 
»un defocus. 
IB. J. T. finch, dau not shown. 

iJ^J"^:.^^^' of the Ttmf johi, infm 
Wosjum (1976 pp. i„S: L ICeU^berjer. K 
naner. M Wwtz. Utuvnknscopy 9, 139 (1982)- L 
S^^re and 0. ^ DeRosier. U 1« 

l i/ '1^3): K V. Nermut 

mott. S. WiUens, S. FuUer. 0. Thomoson. I Biot 
^ifU;!''' 0998): L Bart«lU^?2^Bo/u: 

ivr^Tir*?- *• Acad. So* 

95. 7299 (1998), 
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The response of mammalian fibrobiaste to 
scnjm has been used as a model for studying 
growth control and cell cycle progression ( J). 
Nomial human fibroblasts requiie growth 
factors for proliferation in culmre; these 
growth factors arc usually provided by fetal 
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bovine serum (FBS). In the absence of 
growth factors, fibroblasts enter a nondivid- 
ing state, termed G^, characicrized by low 



metabohc activity. Addition of FBS or purj. 
fied growth factors induces proliferation of 
the fibroblasts; the changes in gene cxpres- 
sion that accompany this proliferative re- 
sponse have been the subject of many studies, 
and the responses of dozens of genes to le- 
rum have been charaaerized 

We took a fresh look at the response of 
human fibroblasts to serum, using cDNA mi- 
croanays rtrprcsenting about 8600 disHnct hu- 
man genes to observe the tcmpoial program of 
transcnption that undcriies this rtsponse Pri. 
mar>' cultured fibroblasts from human neonatal 
foreskin were induced to enter a quiescent state 
by serum deprivation for 48 hours and then 

T^""^^ °^ "^"m containing 

10% FBS (2). DNA micioanay hybritoS^ 

"Tk^^ '° ""^^^ ^ ^""P«^ changes in 
mRNA levels of 8613 human genes (3) at 12 
limes, ranging from 15 min to 24 hours after 
scrum stimulation. The cDNA made from pu- 
nfiwi mRNA from each sample was labd^ 
with the fluorescent dye Cy5 and mixed with a 
common reference probe consisting of cDNA 
made from purified mRNA from the quiescent 



* • » • • 



Rg. 1. The same sertion of 
the microarray is shown 
for three independent hy- 
•'"^Mtois comparing RNA 
isolated at the 8-hour time 
point after serum treat- 
ment to RNA from serum- 
deprived cells. Eadi mi- 
cfoanay contained 9996 
*lwT*ents, indudir^ 9804 
fwjman cDNAs. represent- 
ing 8613 different genes. 
mHNA from serum-de- 
pnved cells was used to 

prepare cDNA labeled with — 

««Lrttaneously hybridized tithe miaaSTC^rS^^o^^t^ ^ "^"^ ««« 

mRNAs are more abundant in the serunvd^ive^ 8«es whose 

as greenjpots and genes whose mRN^ a^ ^^e SS^t*!^^,'^^ 

prote. PS.. 2. ^^r.:!^'^:^,^!^^:^^^^^ 
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cuhurt (tiine zero) labeled with a second fluo- 
rescent dye, CyS The cplar images of the 
hybridization irea^ (Figi iy.wm naidt by 
rep re seniin g the Cy3 fKioTesoem 
fftai afUl,tlie fhicneacein hnage as red and 
mdsing.the two color: imagesL- ^ . Vcc ' : - 

; Dhme temporal profiles, of :^ene' expres- 
sion could be seen amonig tfie 8613 genes sur- 



REPORTS 

veyed in this experiment (Fig. 2); many of these 
, genes (about half) were unnamed expressed 
seqii^« tags diymc 
pattems of expression were observed, Ac, order- 
ly chofepgrqihy of the expression progtam: be- 
caaiK: apparent ^tra Hjc rwulis'^^ analyzed 
by arclustciing aiid.di^iay method developed 
in our laboraiory for ana^-zing genome-wide 



^ CD 



F^ ' 2. Ouster image 
showing the 'different - ^ f t \. 
dasses of gerie e)q)rw^ ""''.o.SS S S £ 
sipn pfofitei Tfve.hurvr 
dred seventeen genes 
::w4x»e:mRNA, levels 
dwged in response to 
senm stimuiation were 
; setecied (7). This sub- 
set of genes was clus- 
tered htersrchicBlly into 
groups on the basis of 
the simiiafity of their 
e)q>ressibn' profOes by 
the procedure of Eisen 
et a/, ley The expres- 
sion pattern of each 
gene in this set is dis- 
played here as a hori- ' 
zontal strip. For each 
gene, the ratio of 
mRNA levels in fibro- 
blasts at the indicat- 
ed time after serum 
stimulation ("unsync* 
denotes exportemiaUy 
growing ceUs) to its 
level in the serum-de- 
prived (time zero) fi- 
broblasts is represented 
by a color, accordir^ to 
the color scale at the 
fcxittom. The graphs 
show the average ex- 
pression profiles for the 
genes in the corte- 
spondir^g -duster* (in- 
dicated by the letters A 
to j and color codir^ 
In every case examined, 
^A/hen a ger>e was rep- 
resented by more than 
one an-ay element, the 
multiple representa- 
tions in this set were 
seen to have identical 
or very similar expres- 
sion profies, and the ; 
profiles correspondir^ 
to these independent 
measurements clus- 
tered either adjacent 
or very dose to eadi 
other, pointing to the 
robustness of the dus- 
tering algorithm in 
grouping genes with 
very similar patterns of 
expression 
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gene expression data (6). An cxamy^ of luch 
an analysts, here applied to a subset of 517 
fo«5 whose ^cpressim chn^ wbsnmiaUy 
in response to serum (7). is shown in Fig, 2. 
The entire. detailed datt m unde^^g^ F^^ 
2 'is>vailable as a tab-delimited table (in 
cluster order) at the 5ci>/ice Web site (www. 
sciencenwg,aij/featuTe/daitily^ In 
addition, the cntnt/ iarger data set for the 
. complete set of genes analyzed in this caqw- 
iment can Iw fwind at a Wdb site iMntaij^ 

by our laboratory (genoauswwwJitanford. 
edu/serum) (5). ' -jr,. i 

lOnc , measure of the; reliabili^' of ;^,the 
changes we observed is inherent in' the ex- 
pression profiles of the genes. For mm genes 
whose expression levels changekl, we could 
sec a gradual change over a feu' time points, 
which thus effectively provided independent 
measurements for almost all of the observa- 
tions. An additional check was provided by 
the inclusion of duplicateand, in a few.cases, 
multiple array clenients representing the 
same gene for about 5% of the genes inchxded 
' in this microarray. In addition, three indepen- 
dent hybridizations to different mtcroamys 
wi* mRNA samples from cells harvested 8 
' hours after serum addition showed good cor- 
relation (Fig.. I). As an independent test, we 
measured the expression levels of several 
genes using the TaqMan 5' nuclease fluori- 
genie quantitative polymerase chain reaoion 
(PCR) assay (P). The expression profiles of 
the genes, as measured by these two indepen- 
dent methods* were very similar (Fig. 3) <yc). 

The ffanscriptional response of fibrcA>lastt 
to scrum was extremely rapid. The hnmediate 
response to scrum stimulation was dominated 
by genes that encode transcription . factors 
and other proteins involved in signal trans- 
duction. The mRNAs for several genes {in- 
cluding c-FOS, JUN B, and mitogen-acti- 
vaicd protein (MAP) kinase phosphatase-l 
(MKPI)] were deiectably induced within 
15 min after serum stimulation (Fig. 4, A 
and B). Fifteen of the genes that were 
observed to be induced by scrum encode 
known or suspected, regulators of transcrip- 
tion (Fig. 4B), All but one were immediate- 
early genes— their induction was not inhib- 
ited by cycloheximide {11). This class of 
genes could be distinguished into those 
whose induction was transient (Fig. 2, clus- 
ter E) and those whose mRNA levels re- 
mained induced for much longer (Fig, 2, 
clusters I and J). Some feanires of the 
immediate response appeared to be directed 
at adaptation to the initiating signals. We 
observed a marked induction of mRNA 
encoding MKPI. a dual-specificity phos- 
phatase that modulates the activit>' of the 
ERKI and ERK2 MAP kinases (72). The 
coincidence of the peak of expression of 
genes in cluster. E (Fig. 2) with that of 
MKPI (Fig. 4 A) suggests the possibility 
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An oontimied activiiy of Ihc MAP kinase path- 
wiy is required to mahnain mductian of these 

genes but na of those with sustained cjqffcs^ 
(dusten 1 and J). The cene 

encoding a seco nd 
member of tf>e diial-qiecificfly^MAP icinase 
^I^Btase hmsfy^ known as duai-spedficity 
protein pho^>hatase 6^p)«2; u« induced^I^^ 
at about 4;houn after serum stimulation. Genes 
encoding diverse other proteins with roles in 
signal transduction, ranging fiom cell-surfece 
recepiors 'ffor ocample. the sphingosine I- 
jAbS|ihate receptor (EDG-I), fte vascular cn- ; 
ddthelial growth fecior rccepii^, and the type 11' 
BMP receptor] to regulators of G^dn sig- 
naling (fff^ cxan^ NET1^II5 rho GEF) to 
DNA-binding transcription feciois, were in- 
<hwed scrum (Fig. 4A). ' 

: Thc reproghrniming of the regulatory cir- 
cuits in response to serum involved not only 
indiiaion of transcription facton but also re- 

duccd expression of many transcriptional reg- 
ulators— some of which may play roles in 
maintaining the dells in G^ or in priming 
them to react to wounding (Fig. 4C). Perhaps 
as a consequence of the historical focus on 

serum stimulation of fibro- 
blasts, the set of transcription facton whose 
expression diminished upon serum stimula- 
tion has been less well characterized. 

Gdics known or likely to be involved in 
controllirig and mediating the proliferative re- 
sponsc showed distinctive patterns of regula- 
tion. Several genes whose produas inhibit pro- 
gtession of the cell-division cycle, such as p27 

Kipl, p57 Kip2, and pl8. were o^iresscd in the 
quicscerii fibroblasts and down-regulated be- 
fore the onset of cell division. The nadir in the 
niRN A levels for these genes occurred between 
6 and 12 hours after serum stimulation (Fig, 
5A), coinddem with the passage of Ae fibro- 
blasts through G,. TTie levels of the transcript 
encoding the WEEI-like protein kinase, which 
IS believed to inhibit mitosis by phosphoiyl- 
ation of Cdc2, diminished between 4 and 8'to 
12 hours after serum addition (Fig. 5A), well 



before the onset of M phase at mind 16 hdini 
raising the'possibility of ari additional role for 
Wcel m an eariier stage of the ceO cvcle W ih 
regulaimg the Gd^ to G, ^tosition:' Swril 
genes mduced in the first few hours' after sertim 
sttmutoion. such as the helii'-liop^lix W- 
tetr^ [D2 ^ ID3 and EST AA0I63O5: a g^ 
with homology to G.-S cyclins; are ciididatei 
tor roles in promoting the exit from Gq. 

Genes involved in mediating progression 
.plough the cell cyclc^eie characterized by a 
distinctive pattern of expression (Fig. ^, clus- 
^ D), refiecting the coincidence of their 
expression with the rwntry of the stimulated 
fibroWastt into the 3 tell-di vision cvcle. The 
stimulated fibroblasts replicated their DNA 
about. J 6: hours aftw serum treatment. This 
timing was reflected by thi induction of 
mKNA encoding both subunits of ribonucle- 
otide reductase and.PCNA. the processivity 

factor for DNA polymerase, cpsilon and delta 
Cyclin A, Cyclin Bj; Cdc2. and CDC28 ki- 
nase, regulators , of passage thfough the S 
phase and the transition from G, to M phase 
were induced at about 16 to 20 hours after* 
^™n addition. The kinase in:ihe Cyclin 
e i-CDK pair needs to be activated by phos- 
phorylation. The gene encoding Cyclin^ie. 
pendent kinase 7 (CDK7; a homolog of A-^. 
nopus MOI5 cdk-activating kinase) was in- 
duced m parallel, with the Gdc2 and Cdc28 
Sr ^jr^- S"fit>csting a poiemial role 
torCDK? m mediating M phase. DNA topo- 
isomerase II a, required for chromosome seg- 
regation at mitosis; Mad2. a component of 
the spindle checkpoint thai prevents comple- 
tion of mirosis (anaphase) if chromosomes 
are not attached to the spindle; and the kinet- 
ochore protein CENP-F all showed a similar 
expression profile. 

hours after the scrum stimulus, one of 
the most striking fearurcs of the unfolding tran- 
scnpnonal program v^-as the appearance of nu- 
mCTOus genes with known roles in processes 
relevant to the physiology , of wound healing 




100.00 ,- 



"1, *^ piognin dioode pmlucis ih« 
cto Assolnoon. and «rn«jdeln,g and jS 

^promote .towtaxis and activitiM bf «o. 
^hfe(for otainple. C0X2) and «^hment 
^, mzyjat,m of monocytes, and maoo. 

for example, inierleukin^. (i^)] 
lymphocyte (for exkmpJe; JCAM I). 

Sg£«*?*. .anrfanugOHsJeciHc 

recruit-- 

ng the phagoc>T.c cells that will be required 

Jwjvound: (ly) promote angiogenesi, aid 
"^r'»f"''>nation, (for: example. VEGF). 
through newly fpm,ine iissue:.-(v) p^m^ 
migration aijd proliferation of fibroblaas (fw 
•*»^P'\CTCF) ai,d their differentiatSSS. 
myofibroblasts (for example. Vimentin); atS 
(vij promote migration and prolifeiaiicMi of 

kOTimocyte^ leading to Veepithelialiiation. 
ofthe wo«,d (fprexample. FGF7). and pro! 
mow prphferation of melanocytei. peij«i, 
contnbutmg to. wound hyperpiimema^ 
(for example, FGF2). 

^ Coordinated regulation of groups of«enes 
whose products act ai different sieps in a' 

example. Fnnn. a prohormone-processinB 
protease required for one of the processiS 
steps in the generation of active endothdiiL 
was induced in parallel with induction of the 
fpl- '"l^iia^'iP^OT of endothelin-l 

|-ALLA/CDIO. a membrane meialloproteaie 
that degrades endothelin-l and. other peptide 
mediators of acute inHammation, .iwTre. 



MtcTMrrayt 



10.00 • — 



1.00 



0.10 



Rg. 3. Ir^dependent verifotion of mtcroarrav auantrbitiftn d-j^- 
teveb of the bx(i£«ed g«»s (Man. n^a/nS^^fSHc'lS^^ 
were measured wKh the TaqMan 5' nudeatt OuJSl!^ receptor) 
assay (9) (left) in .he sam.%Js tte ^e^TLS^""!:! 
micro^ray hybridi«ion. (STo,^ tS1^^S?^5S^'* 



Time 




normalized to mRNA concentrations and plotted rH*tK«. to .k. . ^ 
time zero, so that the results could bico^^ll^^ « 
m-croarray hybridizations. In general quaSEt^th^^^ 
gave very similar results (70) ^"anwiatKm with the two methods 
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Aiced; A iccohd exmspleis provided fcry>a set 
of fivc^encs,mvolvcd,ni the biosynthesis of 
cimlcsierol (Fig. 51), Tlic liiRNAs cto^ 
eadi pitfiae enzymes showdshiuply dimin- 
ish^ CKiweMpntwgimm^^^ hours aftCT 
sewrstinnilaticm pfn 
planatiori for thcrcooWira 
tion of the cholesierol bipsynthciic, pathW^ 
is that icniin provides cholestero( to fibro- 
blasts through Jow-dcnsjiiy .JipoproteinsV 
;»*ereas. in' tlw. abswicc of the choJcsterol 
pipvjdcd ify scnim, cndogCTious cholesterol 
biosynthesis in fibroblasts is rwjuirwl. , • " 
J .Many, of the piwqusly studied gcnes^^^^^ 
c^scn^cd lo.bc regulated . in this pibgrail' 
haye no recognircd role in any aspect of wound 
healing or. fibroblast proiifgratioa Their identi- 
fication in this sriidy maj^ ihmfdn point to 
Previously imlciiown aspects of thcsc'pTOcsscs. 
A few selected genes in this groi^ are shown iii 
^ig. 5K The stanniocafcin gene, for example 
(Fig. 5H), encodes a secreted protein without a 
cicariy identified fimction in human cells {J 4; 
^5). Its induction m semm-stimulated fibro- 



bl^ suggest, the pbs^^ tha, i^^y |^^^: 
role, m thc^ wound-hcaling, procicss,^ pertiiw' 
seryu^.as a signal in' mediating iJiflanmiktion 
oramgic^encsis. . 

';pne. c,f the mos^ inpoitam 
exploratimi^vas. the discovwy of iw:200 pii- 
vio^ly unknown jeri« whosi .cx|iressidn tii^ ' 
regulated in^ specific temporal pancms durini 
the iwponsc of fibroblasts to scnmi/Fdrcxam: 
Plc.;l,3 of the 40.gcnes in chisicr D (Fig.'2) have 
desCTptiA^. names, that .reflw' thcir^p^^^ 

genes;(69V«) encodr 
Protems- thai play rolb in .ali ■QTie.progi^' 
aoiv paniculaxiy in. D?^A Itplicariori mi ni' 
Ga-M iransiiion. TTjis enrichment f6i cell' 
cyclerrelated genes .siiggins tt^i.soiiie of the 



t.t 



9«fe Similaily. a remaHcable fcwkmof gaei^ 

w«,pH,ped iiito cKwerFci, the 
iheir oqprcssi ptweins ini 

vplyed in intei«lhilar sipuling (Rg; 2i, 
eesft^g thata siiniliir role'ihbiild be obhsid^ 

«ae-in this' cluster; a 
d^wopomonrob: large;fciriion of the gend^' 

soipUlanon »-eit waamed ESTs. ™ 

v>rm.m,ion wijstoWth^ 
a model to . stiK^ the. om^ of the. tnniiW ' 




Fig. 4. "Reprogramming^ of fibroblasu. Expres- 
sion profiles of genes whose function is likely to 
play a role in the reprogramming phase of the 
response are shown with the same representa- 
tion as in Fig. 2. In the cases in which a pene 
was represented by more than one element in 
the microarray, aU measurements are shown. 
The genes were grouped into categories on the 
basis of our knowledge of their most likely role 
Some genes with pieiotropic roles were includ- 
ed in more than one category. 



Fold repression 
>S >6 >4 >2 
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Cholesterol biosynthesis 



Fold induction 
*^ >4 >€ >8 
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to seam, stimulation that appear to SI related tf *^,"f"'P"°""' P??'*" of fibroblasts in respom^ 
fibroblast proliferation are sCwifh the same corntEren!!!: "".'"""^ ""^"^ 
teveU as was used in Figs. 2 and 4. (A) Cell cyde anroS.,^ * ''"? »»T»aipi 
(C) inflamniation. (D) anglogenes s (E) ti^up V^ '^^'"^'^^ 

reepithelialiiation. (H unfdentified role in^oun^^^^^ 'eorganiatioaW 
n«nbe„ in (C) and (C) refer to gen« wtese oro^f.«^^ biosynthesis. The 

monocytes and macrophages («,. Vph^^Vs aTi,^!;::^;^^^?^^^ gj. 
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. iCxpresaon;pn)fDmg cxpcrimcms is that the cx- 

• the processes that actually occur 

.'c and not mbely the tingJe process one iiitteded 
-1 *o oteerw. Scnim, the sohiWe fiaction i>f clot, 
tvted blood, IS iKHinally cncoun^ in 
vivo in the comort of a wdund Jnd£cd.^the 
^' CApiiia aion pTogTam that^wc observed in re- 
sponse to senim suggests that fibrobteis ire ' 
: Progrmimed to intcipret the abrupt exposure to 
serum^not as a general mitogcnici stimulus but 
asa^Jec^ ^ 

^votoL The proliferative response that w oric- 
p mally miended to study appeared to be part of a 
*«Fcr physiolppical response of fibroblasts to a 
'•wound. Other, features of iihc transcriptional 
re^Kmse to senim si^ggest that the fibroblast is 
^ ffli.ncdve paipc^ in a convmation anxmc 
-the diverse cells that wk . together in wound 
n^ir. mtcipreting, amplifying modifying,: and 
ibroadcasimg signals controlling inflaihmation 
' angiOgenesis, and epithelial regrowih durinc 
tnc re^xmsc to an injury. 

We recognize that these in vitro results 
almost certainly represent a distoned and in- 
complete rendering of the normal physioloc- 
1^1 response of a fibroblast to a wound. 
Moreover, only the responses elicited directly 
by exposure of fibroblasts to scnim were 
cxammed. The subsequent signals from other 
cellular participants in the normal wound- 
ncalmg process would cenainlv provoke fiir- 
thCT evolution of the iranscripiional program ^ 
in fibroblasts at the site of a wound, which 
this experiment cannot reveal. Ncvenhcless 
;^ believe -that the picnire that - emerged 
strongly suggcstt a much larger and richer 
role for the fibroblast in the orchestration of 
this important physiological process than had 
previously been su$pccte<L 



Repokts 



Refcrcffictf and Notc^ 

(ATCC Oa 2091) in pmg, B w« iJTin 
P«ri dohcs in Dutbecco'i minimum tucmiil 

Zl Th* ceOs wtft thfn.washwJ-ihfw 

"TO w»i ihff umr mt«um iKkm rfis ««l w 

•KWtd from ««v«nl ptMes oT ccUs *«fvwirt btfor. 
»n-T. stwulrtonc this mKNA mv«r«Wi«!«3 
aMtimr.im» M tr«nc» s«nplc. CdU h»- 

wrts (15 Hwx 30 mirv 1. 2. < 6^ a U 16. 20. «id 24 
^L^^JL^ of «rum. mRNA 

-^Jrom .xponmti^y po^ fibrobUiu (not 
a^eaed to serum itjfv«tion). mRNA was isoUttd 
'-lA ttjeFwTnKk mRNA bcbtion Icit (Im^trogen) 

fnwlium was rcmoird. tnd Iht aUs wt,. JuicWy 



■ iWonng hum* $^u-J!^!^^ 
■ >::>79L At ttJti^ JL^^. wwBoiptf or genes (18, 

Of JWOO.of these UnfCwTH^^ 
-• ii«wte«ll\u^ >«™»<«««y from 1 nuMtr ol 

1- t«?«*fOi Cenetia Tha "lOK s#t- it - 
j "^*n. more recent -ISK set" besWet^ 

«*»clonei472are^ita^^i^^^^ ^ 

tpresertt 8141 dinina d^^T^ 
Cerws. m UniCene These or hurmn 

ram.. . -i-TT^ aones. thus presurnMi tn 

8613 different gwes. wcr* o-TrT^ ^ 

erowrew* »i«fh^L^ ^ ™^ hybndaii^ to the mi- 

«0NA (row this b«*!!l! ^^'^ "^y^-"^^ 

Utkw anrf (r*-. tntervjU tfter serum stimu- 

* """Petition. 
Of eenes Hikf^ " *» « i«rge number 
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tumoursfrom which the cell linef were derived The con^^ 1^ ^ ^he;^. - ; 

tumour, stromal and inflammatory components of th« . ^«- " reflecting the 



Introductibn 

Cell lines derived from human tumours have been extensively used 
as cxpmmental models of neoplastic disease. Although such cell 
lines differ from both normal and cancerous tissue, the inaccessi- 
bibty of human tumours and normal tissue makes it likely that 
such cell lines will continue to be used as experimental models for 
the foreseeable future. The Nauonal Cancer Institute's Develop- 
maital TTierapeutics Program (DTP) has carried out intensive 
studies of 60 cancer cell lines (the NQeO) derived from tumoun 
from a variety of tissues and organs'-*. The DTP has assessed many 
molecular features of the cells related to cancer and chemothera' 
peutic sensitivity, and has measured the sensitivities of these 60 cell 
Imcs to more than 70.000 different chemical compounds, includ- 
ing all common chemotherapeutics (http://dtpj)ci juh.gov) A 
previous analysis of these data revealed a connection between the 
panern of activity of a drug and its method of action. In panicuiar 
there was a tendency for groups of drugs with similar patterns of 
activity to have reJated methods of action**^'. 

We used DNA microanays to survey the variation in abun- 
dance of approximately 8.000 distina human transcripts in these 
60 cell lines. Because of the logical connection between the func- 
tion of a gene and its pattern of expression, the correlation of gene 
expTOsion patterns with the variation in the phenotype of the cell 
can begin the process by which the function of a gene can be 
inferred. Similarly, the panems of expression of known genes can 



reveaJ novel phenotypic aspccti of the ceils and tissues studied*^" 
Here we present an analysis of the observed panems. of gene 
^/.STr"" «lationship to phenotypic properties of the 

WceU ImesTTie accompanying repon" explores the reUtionship 
bpveen the gene expression panems and the drug sensitwity pro- 
aes measured by the DTP. The assessment of gene expression^ 
terns m a multitude of ceU and tissue types, such as the diverse^ 
of cell bnes we «udied here, under diverse conditions in vitro and 
w vivs, should lead to increasingly detailed maps of the human 
g«ie expression program and.provide clues as to the physiological 
roles of uncharacterized genes' The daubase. plus toob for 
analysis and visuaLzation of the data, aie available (hnp://genome- 
www4tanford.edu/nci60and hnp://discover.nci.nih;gov) 

Results 

We studied gene expression in the 60 ceD lines using DNA 
microarrays prepared by robotically sponing 9.703 humin 
cDNAs on glass microscope slides"-". The cDNAs included 
approximately 8.000 different genes: approximately 3.700 repre- 

!'o<^uTu'°"*'': human protein's, an additional 

1.900 had homologues in other organisms and the remaining 
2,400 were identified only by ESTs. Due to ambiguity of the iden- 
...y of the cDNA clones used in these studies, we Jtima^ 
approximatriy 80% of the genes in these experiments were cor- 
recuy identified. The identities of approximately 3.000 cDNAs 
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Kfoa 64 cell ttn^I Th. iT •»« of 1.161 tONA, 

4 of «) CM Ibis. Thi. e«««^^itSi?^iL'fL;r^"* "* P«o' in « 

lion l»tl Moii the 60 ctll iin^ («Wi«7hS -.T!^ "li"'*" •»«»"'' 

Mn« KS62 end MO?. The t«o all «„„ ^2^2^^ «P'«enutioo. of „ct, of ihe „„ 
^iSWedfof thesen.durterinsw.i^Xftt^^riH.ili''"* "T* ""-Pondinsly 

onemible t«u. of orisii. of th.^ B^l^^^J «o ref !,« 

pit. pronm; light bluZ tung; orenoe iIk * »»"«f pur- 
bU* unknown OKUAOR-MS)). w'^^ STi ofti'L'.*):;."'*-' 

re^ .oeffitien, rep«.««ed by the Iw-Trt ?i jJ!!!!" T'" *»* 

p.« of node.. Not. th« the two Wplet.^?^Ljte«ILt!n? ''r «onnenin9 
«fl«ly together ^ we,, well ««lff.f^,«l fcS^ 't'l"" •"<• *«<^) ch»u; 

"-.cing thn ,hi. dunging of all li„« htaoS^TSl,!;!."'" """"^ «" »~«. 

reposentnion of the dn. tMe. with the rowi to^^w^if *•• * 

order. The dendrognm apre-nting t^^^^^l^-^'r'"' «"«) » cM-ter 
ted for .lerlty, bu, i. ...ItoWe (httpi^ge^,;^^ ^'^'^ f"" »••« omH- 

all of thi. able reflet, the '^I^^^^^^"lZ1-^:^y^^- *" 
(column). The colour K.le u»d u, re^^^t:*'""*": 'nd cell line 

3..3d-in«,,r.f.r«,the<h«en,*g.^S^^^'^^;"«» » The ..beB 



I. 1 : 



' . • ".',1 ■ 



cell lines 



Each hybriduation compared CyS^labeUed cDNA reverse ,r»n 
-^^ ^r"^ transcribed ftom a refextncc mRm 

inc celi lines (chosen to maximize diversirv in d^t,« ^ - 
d«er™„ed primarily .ftom nvo-dimS:L' d STb^ 
«rnp«mg cDNA from each cell line wiA a coS.?n 
variauon in gene expression across the 60 cell 71 .!i I. 
inferred from .he variation iJuifno^^i ScS 

rauos across the hybridizations. ^^^^y^/CyS 
To assess the contribution of anefactuaJ iounres of variation in 

MCf7 ceU bnes were each grown in three independent cultures 

«tra«ed from each cuhiire. The variance inle trSlicatTfluT 

£J^enT • minimTv^fen h" 

fluorescence signal was greater than' approximately 0 4% of h! 

m«su,able total «gnal dynamic ran^ abml bacW^d ?„ 
either channel of the hybridization We selecTed r 

tie 1^. ^ "rJ!"" to idenjfy 

or I "^Z^"- P'"-'^*' correlation coefficient 

for the triphates of the set of genes that passed this ou°Jh» r„n 
trol 1^ (6.992 spoU included for the MCT7 "mSeland 676," 
spots for K562) ranged from 0.83 to 0.92 (for Si fnd f 
see http://genome-www.stanford.edu/nci60) ^^^'""^ 

To make the orderly features in the data more apparent, we used 
a hierarchical ch«tering algorithm'"* and a pJ^oZo^Z 



i^ig. Ih we concentrated on those genes that showed th^ 
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or1g7nro7S|'c\",C ''''''' ''^ ^^'^ 

Swom" o?'' ^ cell lines 

"ve organ typ„ w.th few exceptions. Cell lines derived fl«™ 
non-small lung carcinoma and breast lumouA were dul^ K ^ 

(Fig 3fl-W). For example, a cluster of approximately 90 ««« wH 
highly expressed in the melanoma-denvVd linTs (S ic) TWs^« 

DCT S IZZor' '-tomera^T?R .5 

s?:=r-^„-— ^^^^^^ 
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i f*" ? <i«».wprcttion patterw rclmd to 
,^^htr.cell-not phtnbtypc«: appJicd 

two diJiiciiiiofMl hierarchical cluttering to 
-nprcsion data from a Mt of (;831 cONAs 

mtasurcd acraa the C4 can lines; The 6^831 
; . cpNAs were iho» with . a ; minimum . f luercf- 
' i*««««fty «< approximately^ 0.4% 

of the dynamic range'above background in 
the reference channel in cadi of the six 

• hybriAiatiom u»d to establish reprbdudbil- 
hjr. Thfi effectiMHy selected, these spots ^hxi 

, proirided - the most . reliable ,ratio , measure* 
mentt. and. therefore identlftcd a subset . of 
penes useful for exploring patterns comprised ' 
of those whose variatioh in expression across 
the 60 eel) Jines was of moderate rnagnltude. ' 

• A Ouster-ordered^ dau ubie. c ObuWing - ' 
: time of cell linet. Cell lines are given in dister 

order. Values are plotted relative to the mean. 
Doubling ttmcs.greater than.the'mcan are* 
shown in green, those with doubling'tim^ tess ' 
than the mean arc shown in red. tf/'Three 
' related gene duRere that were enriched for 
'<gcnes whose expression level varietion was 
correlated with cell line proliferation rate. 
-Each of the three gene dusien (clustered 
solely on the basis of their expression pat- 
terns) showed enrichmern for sets of genes 
' involved in distinct 'functional categories (for 
example, ribosomal .genes versus < genes . . 
involved in pre-ftNA splicing). e. Gene duster 
in which all charaaerized and sequence-veri* 
. fied cONAs encode genes known to be regu- 
' lated by interferoru. f. Gene clutter enriched ' 
for genes that have been implicated In drug- 
metabolism (indicated by atterislts). A further ^ 
property of the gene clustering evident here 
and in Fig. 2 is the strong tendency for redun- - 
dxnt representations of the same gene to 
clwtier immediately adjacent to one another; 
even within larger groups of genes with very 
similar expression patterns. In addttion to 
iltufirating the reproducibility and consis- 
terKy of the measurements, and . providing 
irtdependent confirmation of meny of our 
measurements, this property also demon- 
strates that these, and probably all; genes 
have nearly unique patterm of variation 
•cross the 60 cell lines. If this were not the ^ 
case, and multiple genes had identical pat- 
terns of variation, we would not expen to be 
•b(e to dtsttnguish. by clustering on the basis 
of expression variation, duplicate copies of 
individual genes from the other genes with 
identical expression patterru. 
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mclanonu). UaxiMVl. the leventh line desienat«i ». »-i 

iKJhtedlittna patient S^Z^^J^L^ ^^"^^ re^nding diff«hc« £ iS^I?^'"^ n»y «flea «r. 

c«?m8 tnuishtiin nuduieiv 



MB«5^'ik>lated ftor^pSSr ^""^ ■■ ^"^in this laxK diW^2^ prolAniting cells (MK167). 

p^ytlut thepHgii, p^^ S^^i^^^'^-'^^S^T^^^ 
, m gene expression patiem is reined to the n«.JZ!Jj^^ • , -encoding' RNA^SDlioB. fcZ- ^ ennched 

tu^ofsobeb^astduj^^^^^ » 

' nay have originated fen, aVnebnomS^ ^'^'^^"'ii*"" Ae.«SSSJ?c^S^ 

, tKat:.the,patientl«dico:*^ 054). supporting the notion Sit ^i^r^T'?®"'"'"^ 

. -TT>e high«r-Je*ri oiganization of tlix^ : ^. '^P^Kd in Relation to - " ^ '^"P'^ 

^secellhW ;.:.: ^^^^^f^^ 

Shared biological prop«Ties of the tissues W, ^ u l « ' ""^ ^ smallcrgenc duster (Fie 2ii) ^li «f .w 

...Igs^ed^^ '^'^«"^;.'-4'Pbe^*^^^^ 

into major .branches that s«arated those t)^r«l!«L ^"^^ of interferon-regidated «en« " ^*'«»n«' 

. charaaeristic of epithelial dST&om ^, ^ »™« apression (data notshtiwnyi^ZT^^T" P*"™ 
more typical of stronuJ cellsXSSf^^ ?''^^ *^^^ 

that is'Sos. strongly ^^t'S'^S'',^^^f;f'^■: ^^.^ ^^^i ^ g^bjS„1r.Sr"« 
carcinomas, six of seven ovarian-derived ceU I^' colon Anoth» cluster. (Fil ; 2e) lomS 

breast cancer linespositivefor,h^SL,!c™^"?J P™''^' Possible i^TirtS^I^ '^'^ 

genes in this clos.rha« been irSf^toX?,^*"'"*"'^ "^^^Jing gluuWcysSSIcTSr .u'*™* ^" 
epitheliiUcdIbidogr".TT.eduste^^Si^J^;'' «We for the rate lim,Sep^/«JS 

that encode putative transcriptional resulatois of «;filri*"'* ■ " , : 

Phogenesi^ a human homologue ofTSSo»W!b^ / *«'«tate interpretation of no« 

^ithelial-e^^sed tumour -^r^l^^^^U ^dthoC" cE Spl^'""'' 
S^hfSl^-S^^^^^ -um-mediateiTdtin— ^^ffi^X^r^^^ 

renal-cell.can:inom,^erivedS ^d ,r?mS- *" n,igh, pr^v^a fram^^f ^ 

noma-derived lines. The characterise^ « „f ««a">mg care.- expression panerns in niTour "terpreting the 

this duster indudrf ninTiS^ulLtf ^^f"**^ two breL cancer b ^"^^^^^^^^ ^ """P'"** 

mal cdl funoions (Fig. 3%. InMeS^'SE^o^SX ™' ' "•'"SS^h^s''^^^ 

desenbed as 'sarcoma-like' in appearance (Hs57«tT ""finally cancers (excluding MDA-MB-435 and MnA u , T 

nosa^oma. and SF539. gliosa^ Spl^^^^T" ™' ''""-"8 h!ghliKa^;S;7th"^J^^^ 

genesJ^JS. Although no single gene was^W^Ti ^ P'"*™ shared between the canr„ ^Z"* **P'**- 

of this duster. Z cell iSe'sW^a" d.t™^^^^^^ "lli^"" "^--d from br:!". 

«^*sion of genes encoding proteins with rolesTs^itE or ^""^ '""^'"^ ^ (KRT8) and Sn ig f Kimo, 

modjficauon oftheextraceUular matrix ffo,-,I» i ""^ well as most of the other eDitheli,!' ». ««?tJn 19 (KRT19). 

tenstic of both epithehal cdls and stromal c^.£C ut,' f "'P" originated in tumour ceUs with f«!^r« . . 

dustered with the CNS and lenal «n^!^l„ ^"^'^^y '""""^l epithdial cells (Fie 5^) P^-/ / °f 

because genes characteristiS i^^o'" ^ «'"'«'^ «^ stromal c ^' .ni&r^' "^^enes char- 

nioreabundantlyrepresentedin'Si^J," SS^^, ^°^^) sZ J S^^^^^^^ 

(i/iGLN), was a feature sharwi hv fk • markers 

Physiological variation reflected stromal-Iike cell hnes Hs57ST >n6vlli9%T.T?:. 

in gene expression patterns «^ expression pattern seen in the tumour t™,i 
A duster diagram of 6.831 genes (Fig. 2) is usefii] for exoloH ^"J '° "^"''^ "mponen of h^rm^S" 

:X;;&t:^o=r^^^^^^ 

dusters that were «rid,ed for .n^^^^T^^^^S Tr^u^"' ^^^^^^^ oft SSf 

tne tumour (this was confirmed by staining with anti- 

9»n 
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v-iU5ncr ^o/ CCS IS) 

FIfl. 3 Gene clusters r«l«td to tisiue ch.r.aerinici in th» »ii i r . 
Me«$ .nd M0A4I). < duller of ae.Si htahiLT^ •"•'•''«>'»»-«*'iv.d line* (6/7) .nd?««?:..T^i', " "P'"""! « . low., level i„ .1, J. " 
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l.j. 



lubset of the NOM cell lines to euiore J^SJ^J^Sl^ ^ from • 
tndicJte gene cluiters (shown in dMaii i« »««pm. Ubto 

fUined with tnti-terMin .ntibodiB. rtB^^^'^/^ Hueimtn 16 
tellol., compontnn of tha tiiiue»^n it!7^ ** «««^t 



/•■If,"' * y, 
peeHlmtidn cluster 




te^ mff J^Wn? B«» with rebttd opresrion pat- 

.,'^<^.icn^^ ^rhm Ac clusters. UnchaxaeS « 
Jkely-,0 encode proteins thai haW >„,« snSo tS^f ^ 

W„ gene p^duas with-whid. they appTto W^I.^ 
StJ prsey^iald^^^ 

' • ^ •-^^'*''*"»ypwide insight'' TheieUtio^ii f 

im^unojobulin antibodies; data no, ihown): Th«f6%'ais- '!^%"^^^''^'°'^^^tis«^ryj^,n^r^^ 

tinaseuofgenes wthcp-varyingejcpression among the sai^DlK J?^ »«> example of linking variation in eene 

(fig. 4, aiTow) appear ,o represent, distinct cell typ5 J,at ^'^S T^Z^::^ """"""tie and divers^ phenotypicv^Sid^r 

distinguished ui breast cancer tissue. A fourth cJu$t« of «J^.7 V ^« P*"?'"!' of, gene expression measured in the NC160 c^n 

.nore highly expressed in all of the cell lines Lt rany K SnT'""^ V '^'"'T^'' ''«'P' ««» distin^^i^ h^ 2^ tSS 

clinical specimens, was enriched for genes present in a/w?/ ^ "'S'"" histologicaUy comola b,«^ 

eration- duster described above Jig S^l^c l^ul^^t ''^'i™"*''- i. is now fSe ,oS^^^^ 

«pres,ionofthesegenesUkelyparalleledthedifferen7eTn„^^ "P'"""*" tumour spedmS^^inb^^^^^ 

erationrate between .herapidiy'cydingcuTnirt^^^^^^^^^ Tome" Jk", '''' P«»iW^^^^« a^J^U^m 

much more slowly dividing cells in tissues. „T . I ' '"'"""^ «"Pl« by s«npSem 

' ' ^•'^*""*''"''""n^""onalmorpholo«csd3h^^ 

Discussion might be able to observe interactions^^ a tW d 

Newly available genomics tools allowed us to explore variation in fi^'n"^'"''''!.'!'"'"','" «»ationships will bSri! 

gene expression on a genomic scale in 60 cdl lin«"Sfa.™ I t ^ *""»'''""^'y*« of gene expression panemVfrom Ln « 
dh*.«.,umourtissues.VAteuseda,simpleclC«d^^^ weU as dissected tumo„rs"H..s.4.^ . P fro" '"taa a, 

tify the prominent features in ,the gene expression iMtterns thai 

appeared to reflea >olecular signatures' of the tissuTf^^ Methods 

which the alls originated. The histological characteristics of the «»« 9.703 human cDNA clon« {ltee.„h C«U ^ 

ceD>es that domi,«ted the dustedng w^e pervasive «ough puL^ " " """^ co.o„« i^SSSir; 

i^?w T"'"! >*en alternative subsS Srunt™) w ' '«P««"-ing?o^ 

SpSnTnet-n^y £ Afi^^oll^Sf^^^^^^^^ '^^^^^^^^ 

have presumably all been shaped by seleaion for resistance to """'"y ' "I"'"" ««iiot obt.ined we .iien,™^ fwr""* 
ho« defences and chemotherapeutics and for rapii S»non ^'"''"^ °^ «f 'SS ° ""^^^ 
.n the tissue culture enviromnent of synthetic grLh^SlZ ""'^'^'"f -.ccing p,„.„, of exp,c^„„ Jsg'^.'^^,;^""" 
bovine serum and a polystyrene substratum But rhl „^ incorrectly identified, and soo ind«i™ ?7 

identifiable faaor >ccZiZ, for varia^S ^^«e^S Sna r"'""'^"'*™'''*'''" '''«<^°^^^^^ 
patterm among these 60 cell lines w„ the identity ofXSue -3 cS:"v? ^ " "'^J'"' ^' ^ Sf« S 
from which each cell line was ostensibly derivrH r u Ji'T ^' f"" °f done* used and iheir nrmi 

ceD lines we examined. naihT^^Z^"^„I°l^' ^^'ZTJ^^^r'''': """" P^^^^.W^hX^ZS::-^^' 



W«,. Neverthde,;. the prominenc^ of re^^'dT^ ^^^^ "! 
.he ceU hnes isolated from glioblastomas and c^nom« m J 
reflect a selection for ihe relative ease of establishmenT"fT3 
hn« expressing stromal characteristics, perhaps combined with 
physiological adaptauon to tissue culture conditions'*-'*'. 
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(100 Ml) WIS purified by cd odusion. eant^f^*^ j 

5XSSC no JT-Th^ PTU »,3Tr «nd rejuspended in 



, of mKNA and refcraiKr pobL tm*^ ' " T .^^ . . 
(2 mM) «id 5% f«d «If .o»m.;ToSS,^^„ "d. gh««ninc 



.7>c ; 



.nalyii, were quicuTfroa^ iS ZL^JT'?!! "^'^'W 
«|.ll,pi«« (-5(^100 mg eich). inJSji*^ jj* «"» 

■ inotuing ,0 .20.000 r.p.m..<«« . Stf iffi vi^"^ 

Ml/.umour homogma,e - docriW ^Ae^, ^ ''T^ ™- 

«.j.u„g mRNA «.„i„g &oB^SKTh!!z^^ir'' ""r" 
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Hg. 5 HiftoteBic fc»tur« of bfe«t cancer biopsie* «n h. r-r^. 

dii.9r.m in Fig. 4 shomng gene ciuiten enriched for oem^ J**'**** °" 9«« «»pre«ion patterns Enlaroement. «f ,k 

tuhored cell linei. A climei includins many oef>ei chl™! • . <»t«erem cell types in the bre.it cancer soeciLm f * 

the oestrogen receptor and tumours^. ^^'es ^.^^'^^V"^^ 'i'^*"*' «-Preued in cell "« (1470 ^d Mc Al dl^!!:"?"'*'*!^ ^'^^^ 

specimens. Expression of these genes in the xurZ^^^^^ eX^Tthr" ^''^ ch^-nl'^^^^^^^^^ '^^^^^ <or 

^^^^^^ 
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0VCAS.3 «id OVGAR-4 <o*Mi«a);,eAW.| (ren»n- PC Vf^LiT; « J 



•« hitp-y/nn.jwnfoid.«du/ioftwi«). Each ipot mi dcfiMrf h. — . 
P««..omng of . grid of cWa cm A, 

locd b«kgn««d «« compuMd lb,.«eh .pot «^ to T^^Tl 
miemny „ . upurr of 40 pad. i„ widU, .ml cmwd«X^ 
cen,^. eluding .D pod. wiAin «v defined .pTNrS^^T 

-KT- "V- - r' "««umi ceil popuhUoni using lulforho- int^i,^!^ r^?""^ ^ "Pplyingi, unifom Kak ikctor to aU 
d,muKBopncJdeMhyme.su«menti«y.TT,edou Si^.K clunnet The nom,Ja,tion teor *2 

w« dcuta^d wing the equitioa- N/No,= ek., ^ iJiZ.^"^^ ^i^ "ojfCyJ/CyS) for" i;'.„»ii« „/.pou Z 



SjSitS doubling .im^^baied on ,ou- 

Wrag tuna for crilj inocukted mte 96-wdl platei it the u^i^. • ' i : 
Uon de«i.i« .nd g„,^ in RPMl ie*0,r,Ll'Z^n^':^ 



rt«r48-h moitauon. uxJ t is 48 h.n» eqwUon was then uied with 
|te?ved k to cJcuIm the doMUing tihte . by setting n/No - Z F^g2 S 
K-rob.«n^ and N value, by" .^gi„g dewitie, 

.... 



.1 1.* 



V) 



Prep««tion and hybridiiatibn of fluorescent labeUed cDNA. For ad. 
coraparaiive array hybridization. labeUed cDNA was nnthesized bv r^JZ 
iramcripUon from test cell mRNA in the presence of CyS^dim-. 
the rrierencenJUM with Cya-dUIT. uang theSuperscrioTll r«Zf . 
Kription kit (Gib«>.BRL). For ead, reve.^ .r.n;;^^''^;^^^^ 
(2 Mg) was miaed with an anchored oligo-dT (d-20T-dlArrI ^'-^^ 

g)in..otalj,.„™eof,5plheatedtclox:f:',fl'S^^ 
To thu^ample. we added an unlabelled nucleotide pool <0« ul- « t!," 
each dAT?. dOT. dCTP. and 15 mM dTTT). eithrc^r ^ 
dUn. (3 Ml; I mM; Amerrftam). ^^-^^T^K^Zl^^^'^ 
HO. pH 8.3 375 mM KQ. .5 mM MgQ,). 0.1 M dA (3^?°.n^ iTf 

^ the UNA wa. d^raded by adding 1 N NaOH ( I J m1) ,„d i«:„batina a, 
70 'C for 10 mm. The mixture wu neutrdiaed by adding of 1 NHa7i \ 
Ml), and the volume brought to 500 Ml with TE (lb mW\^ 
We .dd«J Cotl humin DMA (20 Me Cibco-BRL). al^ p^riK^Si' 
by centHfugationin a <>ntricon-30 micm-concentrati^ (Ami»n)^e 
two separate probe, were combined, brought to a volume of 5^ ML «d 
concentrated again to a volume of les. than 7 ul wT^rf , , 

poly(AM»NA (1 Ml; Sigma) «,d tRNA (10 M^rt g1^b£?11 .dTeS^ 
and .dj»,ed the volume to 9.5 Ml with dS^^'^t'T^.'^^i 
preparation. 20XSSC (2.1 Mfc l^SMNaQ. 150 mM 1^^^ J^. 
.0% SDS (0J5 Ml) w^eadd«l to a total final vdl'Jo^^y^^f!^"^ 

^Lt"'":? ^ '"«' -C incubated a. 37 ^ 

20-30 mm. and pUced on the array under a 22 mmx22 mm glass coverriiD 

Z 'S"" f - '^'"8^" " fo"^! 8 h in a custoKS^S-" 
ber wuh humidity mami.med by a small reservoir of 3>cSSC Ana« w^r 
washed lor submenion and agitation for 2-5 min in 2xSSC with 0 1 % SDS 
foBowed by 1«SC and then 0.1XSSC The array, were ".pun diy- 
.nhigation for 2 min in a slide-rack in a Beckman 05-6 taWetopTenXe 
m Microplus urrien at 650 r.pjn. for 2 min. «.">irinige 



applied to sdeci subseu of genes fn)m.:ihe 9^3 £bN^d™«,!^JT 
arrays Before du«ering and disrim. the JoiiiAh^'^K ** 

.U ratio. mea«,ed for that gene. TTie e»trii«SS£«^rT 

we applied • hierarchical clustering algorithm seiMr.t»lC i» .k-Vi i- 
and gene. u«ng the tarson cortdati^n ^^^^^T^^,!^ 
larity and average linkage Amxm,,'^.»-^^^t2^JT 
.wo dendrogf.,^. ,.ree,).-one f^e cen'li^^ a^rfoMhrr* 
whid, very simaar elemenu are connect^nTAort brenAei a„d^." 



>^ ,u«Uution and data proceasing. Following hybridiaaUon. arrays 
were scanned using a la»r-scanning microscope (ref. 17- htmJ cmoJ^ 
...nford.edu/pbrown).Separateim,ge.w«eacqC[,edf„Q^^^^^ 
earned out data reduction with the program ScanAlyte (M B.L. a^^Uaw! 
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»!.. xnea U. rt «/. A gtne expression dcutoise for th» ^ML^.a» " " ' 

-e»neer. Mature GerSt. 24; 23fr-i44 OOoS^ Ph-rmKcrfogy of^ 
12. Khaa, L M aAtGene, exprcnion prof ilitw of alMBft>r *kW*«j L 

Oer. $.D ' Zhou. W!iii«fm ^ « b * .^i^llir*/ 



gene. Ieth««2)9i«m * l«fv.« tlStoi)"b; f^^^Ilr ^ »uppf«ibf ^ 
•PrtheM celis during DrosophiU^^Joc^Sr^^^ chan^eof: 
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